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1.0  INTRODUCTION 


The  purpose  of  this  report  is  to  present  a concept  along  with 
supporting  theory  for  a magnetic-guideway/guided-projectile  system 
that  has  been  devised  for  the  AEDC/VKF  Aeroballistic  Range  G 
(Ref.  1).  The  concept  has  evolved  out  of  a somewhat  broader  effort 
to  provide  alternate,  or  "backup"  technology  to  that  presently  under 
development  wherein  the  launch  sabot  is  retained  during  flight 
and  serves  as  a shoe  for  mechanical  containment  of  the  projectile 
between  guide  rails.  The  objective  of  either  approach  is  to  pro- 
vide a guided  hypervelocity  vehicle  that  will  better  enable  study 
of  the  ablative  and  erosive  environments  that  high-speed  missiles 
and  space  vehicles  often  must  survive.  To  adequately  perform 
such  studies  in  aeroballistic  ranges  it  is  required  that  the  flightpath 
be  predetermined  to  enhance  data  from  wayside  photography,  etc. , 
and  to  provide  for  flight  into  deceleration  devices  where  the  projec- 
tile can  be  retrieved  for  postflight  measurements.  Potential  advan- 
tages to  be  gained  by  use  of  the  magnetic  method  include: 

1.  Elimination  of  sliding  surface  contact  and  associated 
wear,  erosion,  and  gouging  of  the  shoe  and  guideway, 

2.  Elimination  or  reduction  of  problems  associated  with 
shock  interaction  between  projectile  and  guideway, 

3.  Less  critical  alignment  and  smoothness  requirements 
for  the  guide  way,  and 

4.  Periods  of  essentially  free  flight  with  no  track  or 
other  obstructions  to  interfere  with  wayside  photo- 
graphy. 

As  will  be  seen,  the  concept  is  closely  related  to  presently 
accepted  notions  for  repulsive  magnetic  levitation  and  guidance 
of  ground  transport  vehicles  and  rocket  sleds  (Refs.  2,  3,  4,  and 
5).  In  these  cases,  the  dominant  scheme  consists  of  attaching 
high -intensity,  superconducting  magnets  to  the  vehicle  and  pro- 
viding roadways / guideways  of  aluminum  or  other  normally  con- 
ductive materials.  By  virture  of  vehicular  motion,  the  fields  of 
the  magnets  induce  and  react  against  eddy  currents  in  the  roadway/ 
guideway  with  the  resulting  forces  experienced  by  the  magnets 
being  used  for  levitation  and  guidance. 
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The  present  concept  is  similarly  based  on  motion-induced  electro- 
dynamic interaction  but  with  reversal  in  location  of  the  key  elements; 
i.  e.  , the  guideway  provides  a range-length  inducing  field  and  conduc- 
tive material  is  integrated  into  the  projectile  with  the  forces  it  experi- 
ences being  utilized  for  guidance.  Resort  to  this  approach  has  been 
necessitated  in  the  present  situation  because  of  launch  considerations: 
i,  e. , projectile  velocities  of  interest  for  abalation  and  erosion  tests 
are  4,  000  to  6,  000  m/ sec  and  to  achieve  such  in  the  VKF  Range  G a 
projectile  must  withstand  accelerations  up  to  2 x 10 5 g's  in  a staged 
powder -hydrogen  gun.  As  reported  in  Ref.  6,  superconducting  mag- 
nets are  not  strong  enough  structurally  to  withstand  this.  The 
alternative  of  permanent  magnets  would  not  provide  a field  of  suffi- 
cient intensity  and  normal  electromagnets  would  require  a power 
supply  which  makes  them  unsuited  for  incorporation  into  the  projec- 
tile. An  electromagnetic  guideway  with  only  normally  conductive 
material  in  the  projectile  as  proposed  is  technically  feasible,  how- 
ever, and  its  practicability  for  the  case  of  Range  G is  strongly 
suggested  by  the  nearby  Tunnel  F power  supply  (Ref.  7).  From  this 
flywheel /inductor  stored  energy  source,  power  pulses  can  be  made 
available  at  currents  up  to  10°  amp  for  energization  of  the  guideway 
during  a projectile's  flight. 


2.0  SYSTEM  DESCRIPTION 


As  previously  stated,  the  system  concept  is  based  on  motion- 
induced  electrodynamic  interaction  between  conductive  material 
that  is  to  be  deployed  in  projectiles  and  a stationary  electromag- 
netic guideway  that  is  to  be  provided  around  the  desired  flightpath. 
More  specifically,  the  conductive  material  is  to  be  in  the  form  of  a 
shell  around  the  projectile's  forebody  and  the  guideway  is  to  be  such 
that  as  a projectile  moves  off  course  it  will  be  exposed  to  increasing 
magnetic  field  intensity  and  as  it  moves  downrange  it  will  be  exposed 
to  spacewise  modulations  in  this  intensity.  Because  of  high-frequency 
shielding  effects,  the  forebody  will  appear  diamagnetic  to  the  modu- 
lations and  experience  forces  such  that  its  preferred  position  pf 
equilibrium  will  be  the  guideway  axis  where  the  field  is  zero.  These 
"positional"  forces  are  to  modify  the  projectile's  angle  of  attack  so 
as  to  enlist  aerodynamic  forces  that  are  capable  of  maintaining 
flight  along  the  guideway  axis.  Flight  stabilization  is  achievable  by 
virtue  of  the  fact  that  lateral  motion  of  the  conductive  shell  will  be 
through  field  nonuniformities  such  that  attendant  eddy  current  losses 
will  tend  to  dampen  or  retard  such  motion. 
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As  initially  envisioned,  and  as  can  be  envisioned  for  present 
analytical  purposes,  the  guideway  field  is  to  be  generated  by  two 
superimposed  sets  of  four  currents  that  are  symmetrically  oriented 
about  the  range  centerline  or  z axis  as  shown  in  Fig.  1.  As  indicated, 
one  set  is  composed  of  range -length  line  currents  of  magnitude  I;q 
while  the  other  is  composed  of  transposed  line  current  segments  of 
magnitude  Ip  and  length  a.  Also,  as  indicated,  there  are  two  go  and 
two  return  currents  in  each  set  with  those  of  like  direction  being 
diagonally  opposite  each  other.  The  magnetic  field  of  the  Ip  currents 
is  independent  of  distance  downrange  as  indicated  in  Fig.  2 whereas 
that  of  the  Ip  currents  is  of  opposite  polarity  in  adjoining  segments 
as  indicated  in  Fig.  3. 

In  both  cases,  as  should  be  expected  from  the  Biot-Savart  law, 
a so-called  null-fLux  field  is  produced  whose  intensity  is  strong  near 
the  currents  and  decreases  to  a zero  or  null  value  at  their  axis  of 
symmetry.  Anticipating  that  the  Ip  field  will  be  stronger  than  the  Ip 
field,  the  combined  results  of  Ip  and  Ip  will  be  a null -flux  field  whose 
intensity  is  modulated  with  respect  to  z.  Later  analysis,  however, 
will  consider  the  two  independently  with  the  Ip  field  being  treated  as 
the  means  for  generation  of  lateral  damping  forces  and  Ip  being 
treated  as  the  means  for  generation  of  lateral  positional  forces. 

Insofar  as  the  projectile  itself  is  concerned,  the  shape  and  struc- 
ture shown  in  Fig.  4 have  been  adopted  for  investigation  of  theoretical 
system  feasibility  herein.  Whereas,  as  will  be  discussed,  prospects 
exist  for  improvements  from  other  viewpoints,  this  configuration 
yields  an  agreeable  analytical  model.  This  is  primarily  because  the 
shell  on  the  forebody  is  in  the  form  of  a cylinder  and  thus  lends  itself 
well  to  electrodynamic  analysis.  It  is  also  to  be  appreciated  that  the 
projectile's  overall  outer  surface  is  in  the  form  of  a flared  cylinder 
(with  optional  nose)  which  lends  itself  to  aerodynamic  stability  and, 
when  in  aerodynamic  flight,  develops  side  forces  as  a function  of 
angle  of  attack.  Moreover,  by  being  a relatively  standard  aero- 
dynamic shape,  much  data  is  available  from  which  to  predict  flight 
behavior.  The  fact  is  that,  in  the  past,  use  of  projectiles  with  this 
general  overall  shape  has  been  common  in  ballistic  range  work. 
Through  use  of  strippable  sabots  this  shape  adapts  well  to  present 
launch  techniques  and,  as  is  seen,  is  quite  flexible  insofar  as  the 
type  of  nose  test  article  it  can  carry. 
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3.0  ANALYTICAL  DEVELOPMENT 


3.1  GUIDEWAY  AND  SHELL  INTERACTIONS 

For  the  purpose  of  developing  theoretical  aspects  of  the  system 
electrodynamics  to  the  point  where  guideway  and  shell  parameters 
can  be  chosen,  an  idealized,  two-dimensional  mathematical  model 
will  be  utilized;  i.  e.  , the  cylinder  is  of  infinite  length  in  the 
z- direction,  conditions  do  not  change  with  respect  to  z,  etc’  The 
concept  being  pursued  is  compatible  with  such  a model  if  the  following 
conditions  can  be  met: 

1.  All  guideway  current  elements  that  lie  in  the  z direc- 
tion are  long  compared  to  the  radial  dimension  of  the 
guideway.  Despite  the  use  of  long  current  elements, 
the  needed  high-frequency  effects  will  be  created  if 
projectile  velocity  is  sufficient. 

2.  The  shell  axis  remains  approximately  parallel  to  the 
z axis  at  all  times.  This  is  possible  if  in  actual 
flight  only  slight  deviations  from  such  conditions  is 
required  to  develop  needed  aerodynamic  forces  for 
maintenance  of  course. 

3.  All  eddy  currents  in  the  cylindrical  shell  (except  for 
those  in  the  end  surfaces  which  are  to  be  ignored) 
flow  approximately  in  streamlines  parallel  to  the 

z axis.  This  will  be  reasonably  true  when: 

(a)  the  vector  potential  of  their  inducing  field, 
i.e.  , the  field  of  the  guideway  currents,  lies  in 
this  direction  (see  condition  (1)  above); 

(b)  the  shell  has  a large  length-to-diameter 
ratio;  and 

(c)  the  end  walls  are  made  thick  for  decreased 
opposition  to  current  flow  across  them. 

Later  application  of  two-dimensional  analytical  results  to 
design  of  an  actual  system  will  demonstrate  that  compliance  with 
these  conditions  is  quite  feasible. 

3.1.1  Conditions  for  Diamagnetic  and  Damping  Behavior 

The  diamagnetic  and  damping  effects  that  are  to  be  employed 
are  both  dependent  upon  eddy  currents  being  induced  in  the  conductive 
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shell  as  a result  of  its  motion  through  the  guideway  field.  Diamag- 
netic behavior  is  synonymous  with  the  shielding  of  the  shell’s  interior 
which  is  achievable  via  such  currents  and  is  related  to  the  inductive 
(or  reactive)  energy  in  their  magnetic  fields;  i.  e.  , the  induced  eddy 
current  field  required  to  completely  nullify  or  cancel  an  inducing 
field  within  a body  is  identical  to  that  field  which  would  be  required  if 
cancellation  were  by  induced  diamagnetic  polarizations.  In  both  cases, 
positional  forces  exist  on  the  body  that  are  proportional  to  the  change 
in  induced  field  energy  with  respect  to  position.  For  the  case  of  the 
shell,  damping,  or  retardation  of  motion,  is  related  to  the  rate  at 
which  such  energy  is  dissipated  in  the  shell  as  joulean  heat. 

Quite  obviously,  only  positional  and  damping  forces  that  act  in 
the  lateral  direction  are  useful  for  guidance  purposes.  The  fact  is, 
however,  that  to  achieve  radial  positional  forces  by  the  scheme  being 
pursued  it  is  not  possible  to  avoid  altogether  the  development  of  both 
positional  and  damping  forces  that  act  in  the  axial  direction.  Nor  is 
it  possible  to  develop  the  desired  radial  damping  forces  without 
creating  some  level  of  undesired  positional  forces.  Fortunately,  the 
axial  positional  force  acts  first  in  one  direction  and  then  the  other  as 
the  projectile  passes  through  the  transposed  guideway  so  as  to  have 
insignificant  net  effect.  Moreover,  as  will  be  seen,  it  is  possible  to 
choose  guideway  and  shell  parameters  such  that  desired  radial  forces 
will  dominate  those  that  are  not  desired. 

With  its  radius  and  length  being  a and  h,  respectively,  visualize 
for  the  present  that  the  cylindrical  shell  lies  concentric  with  the 
z axis  of  a cylindrical  coordinate  system  where  points  in  space  are 
defined  by  r,  9,  and  z as  shown  in  Fig.  5.  Visualize  further  that 
from  external  sources  there  exists  an  inducing  field  whose  magnetic 
vector  potential,  "Sq,  is  parallel  to  the  axis  and  whose  magnitude  is 
varying  with  respect  to  time.  Because  the  shell  is  conductive,  eddy 
currents  will  flow  whose  magnetic  vector  potential  can  be  repre- 
sented by  A'.  The  total  induced  electric  field,  E,  is  then  composed 
of  that  due  to  the  original  inducing  field  plus  the  reactive  contribution 
of  the  eddy  currents.  From  the  laws  of  induction  by  Faraday  and 
Lenz  it  can  be  expressed  as 


Total  induced  E = — (A.  + A') 


(1) 


where  t is  time.  Current  must  then  flow  along  the  surface  in  accord- 
ance with  Ohm's  law  which  can  be  written  as 


E = A 


(2) 
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where  i is  the  current  density  along  the  surface  and  c is  the  area 
resistivity.  Finally,  in  absence  of  displacement  currents  the  vector 
potential  and  current  density  must  be  interrelated  through  Ampere's 
law  for  all  points  in  space.  With  ju0  being  the  permeability  of  free 
space  and  the  shell,  this  relationship  can  be  expressed  as 

V*A'  = Ki  (3a) 

For  points  in  space  that  are  not  part  of  the  surface  there  can  be  no 
current  flow  so  that  for  these  points  the  relationship  reduces  to 

V2A'  = 0 (3b) 

Equations  (3a)  and  (3b)  will  be  recognized  as  those  of  Poisson  and 
Laplace,  respectively. 

Recall  now  that  in  Section  3. 0 the  ratio  of  the  shell's  length  to 
diameter  has  been  postulated  as  being  large  so  that  two-dimensional 
solutions  for  the  conditions  of  Eqs.  (1),  (2),  and  (3)  are  adequate. 
Fortunately,  the  two-dimensional  situation  has  been  treated  at  some 
length  by  Smythe  for  the  case  of  a timewise  - sinusoidal  inducing  field 
(Ref.  8,  p.  400).  The  treatment  expands  the  vector  potential  of  both 
the  inducing  field  and  an  assumed  eddy  current  field  into  a series  of 
circular  harmonics.  For  polar  positions  of  r < a that  series  obtained 
for  the  inducing  field  is 


- A “ 

A.(r,0,t)  = k 1 CnrDa  n cos  (n0  + §n)  cos  cot  (4) 

where  £ is  a unit  vector  in  the  z direction,  n is  the  degree  of  the 
harmonic,  Cn  is  the  coefficient  of  the  n^1  harmonic,  fin  is  the  phase 
angle  for  the  harmonic,  and  u is  the  angular  frequency.  Pro- 
ceeding to  match  boundary  conditions  under  the  assumption  that  each 
harmonic  decays  independently  and  exponentially,  several  pertinent 
solutions  that  were  obtained  in  the  treatment  are  as  follows: 


1.  Area  density  of  eddy  currents 


4a,0,t) 


k2 

£ n C cos  ( nd  + <5 J cos  f cos  (tot  — i ) 

V-a*  n=l  " on  n 


(5) 
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where  en  represents  a phase  lag  between  Aj^  and  A'  and  is 
given  by 


2.  Vector  potential  of  the  eddy  currents  at  r < a 

- >.  00 

A (rT0,t)  = -k  £ C rna'n  cos  t cos  (n6>  + S ) cos  (cut  - c ) /7\ 

tl=  1 ° 11  Q Q \ 1 / 

3.  Ratio  for  the  harmonic  of  the  vector  potential  of  the  eddy 
currents  to  that  of  the  inducing  field  at  r = a 

j 

| A'(r.0,t)l 

—= ! - cos  c , . 

| A.(r,0,i)|  «*  (8) 

4.  Total  vector  potential  at  r < a 


Aj(r,0,t)  + A'(r,0,t)  = -k  I Cn  rna‘n  sin  en  cos(n#  + §n)  sin  (cut  - fn) 

5.  Shielding  effect  of  shell;  i.  e.  , ratio  of  term  of  the  vector 
potential  inside  the  shell  as  opposed  to  that  with  the  shell 
removed. 

"n  = |A5n(r.0.L)|  = 5I"'n  (10) 

6.  Time  averaged  real  power  per  unit  length;  i.  e.  , time  aver- 
aged  rate  of  energy  dissipated  per  unit  shell  length 

Preal  = 2rr<:  a>'1  S 1,2  Cn  cos2  fn 

n=1  (11) 

Now,  in  addition  to  real  power,  it  is  instructive  to  think  in  terms 
of  reactive  (imaginary)  and  total  apparent  power.  From  Eqs.  (1)  and 
(3)  it  is  clear  that  the  portion  of  induced  electic  field  attributable  to 
reactance  of  the  eddy  currents  is 


d A ^ ( a, 6 , 0 

Reactive  E z (a,#, t)  = — 


(12) 
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A time  average  reactive  (or  imaginary)  power  per  unit  shell  length 
can  then  be  defined  as 


P. 

iraag 


2 n <?[A  ' (a,0,t)] 

37 


iz  {a,0,t)<i0 


Combining  Eqs.  (5),  (7),  (12),  and  (13)  yields 


(13) 


P: 


= — 1 n Cj 

Vo  n=l  r 


imag 


COS“  f 


Combining  Eqs.  (11)  and  (14) 


p 

p 


n real 


n imag 


2g  n 


“ tan  *n 


(14) 


(15) 


The  total  apparent  power  per  unit  shell  length  for  each  harmonic 
would  then  be 


n total 


n imag 


(16) 


Finally,  combining  Eqs.  (15)  and  (16) 


and 


n real 
a Lota) 


2 2 2 
PcP  a 


V 1 + col‘  ^ 

T n 


= sin  f , = S (see  Eq.  O0» 


P 

n imag 

_ 

□ total 


- yf  1 + lan2<n 

2 2 2 
/l0a  O2 


(17) 


= cos  (n  = I A'|/lAj|  (see  Eq.  (S» 


(18) 


The  ratios  Pn  rgal^n  total'  ^n  imag^^n  total'  I ^ n \/\ ^ni I anc^ 
are  plotted  in  Fig.  6 as  functions  of  the  parameter  ^0au(2  n <:  )-1, 
i.  e. , - cot  en.  It  should  be  noted  that  as  this  correlation  parameter 
approaches  zero,  power  becomes  totally  real;  i.  e.  , magnetic  energy 
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of  the  eddy  currents  is  dissipated  as  joulean  heat  at  the  same  rate  it 
is  generated.  At  this  extreme  it  can  also  be  seen  that  the  vector 
potential  of  the  eddy  currents  becomes  small  as  compared  to  that  of 
the  inducing  field.  Unfortunately,  this  implies  correspondingly  low 
values  of  eddy  currents  so  that  to  dissipate  significant  energy  for 
damping  of  motion  the  inducing  field  will  need  to  be  strong. 

It  should  also  be  noted  from  Fig.  6 that  as  the  correlation 
parameter  becomes  large  the  power  becomes  essentially  totally  re- 
active and  there  is  almost  perfect  shielding  of  the  shell's  interior. 

It  should  be  understood  that  such  reactive  power  represents  energy 
transfer  to  and  from  the  magnetic  field  of  the  eddy  currents.  From 
Lenz'  Law,  the  flux  of  this  field  is  in  direct  opposition  to  and  tends 
to  cancel  that  of  its  inducing  field  within  the  shell,  hence  the  shield- 
ing or  artificially  diamagnetic  effect. 

Fortunately,  again  from  Fig.  6,  it  should  be  noted  that  as  the 
correlation  parameter  is  increased  by  a factor  of  approximately  2 5, 
i.  e. , from  a value  of  0.  2 to  5,  the  eddy  current  power  moves  from 
approximately  98-percent  real"  to  98-percent  reactive.  * Note  also 
that  for  a given  n all  quantities  in  the  correlation  parameter  except 
angular  frequency  are  fixed  when  a cylindrical  shell  is  specified. 
Advantage  is  to  be  taken  of  these  facts  in  the  guideway  design.  On  the 
one  hand,  the  interval  between  guideway  transpositions  as  shown  in 
Fig.  1 will  be  established  such  that  high-frequency  eddy  currents  will 
be  induced  that  are  >98 -percent  reactive,  and  shielding  or  diamagnetic 
effects  are  large.  In  this  way,  needed  lateral  positional  forces  can  be 
generated  with  minimum  retardation  of  axial  motion.  On  the  other 
hand,  the  flightpath  will  be  tailored  so  that  excursions  across  the  guide- 
way centerline  will  result  in  lower  frequency  eddy  currents  that  are 
> 9 8 -percent  real  where  damping  effects  dominate.  This  will  allow 
needed  lateral  damping  forces  to  be  generated  with  minimal  extraneous 
positional  forces. 

3.1.2  Diamagnetic  Force  Relationships 

Through  use  of  the  principle  of  superposition,  it  is  possible  to 
relate  the  conductive  shell  and  the  guideway  field  modulations  produced 
by  the  Ip  current  segments  of  Fig.  1 independently  of  the  basic  field 
generated  by  the  Ip  currents,  and  vice  versa.  Therefore,  considering 
only  effects  of  the  transposed  field,  one  can  visualize  that  as  high- 
speed projectiles  move  downrange  they  are  exposed  to  a null -flux 


Intended  to  imply  magnitudes  that  are  98  percent  of  the  apparent 
power. 
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field  of  first  one  polarity,  then  another,  etc.  So  that  the  shell  will 
appear  diamagnetic,  the  frequency  components  of  this  effectively 
alternating  field  are  to  be  such  that  the  interior  of  the  shell,  if 
appropriately  designed,  will  be  shielded  (see  Fig.  6). 

Now,  as  will  be  seen,  the  length,  a,  of  the  line  current  seg- 
ments is  to  be  quite  large  as  compared  to  their  individual  distance  b 
from  the  range  centerline  (see  Figs.  1 and  3).  Therefore,  during 
flight,  the  conductive  shell  spends  most  of  its  time  between  trans- 
positions where  a null -flux  field  exists  whose  strength  is  for  all 
practical  purposes  equivalent  to  that  which  would  exist  if  the  current 
segments  were  of  infinite  length.  This  can  be  seen  through  applica- 
tion of  the  Biot-Savart  law  to  a simple  set  of  segments  of  alternate 
polarity  which  lie  on  a straight  line. 

Consider  the  situation  shown  in  Fig.  7 where  such  a set  lies  on 
a z'  axis  parallel  to  the  z axis.  Let  r ' be  the  normal  distance 
from  the  z ' axis  to  a point  of  interest  and  n stand  for  the  n^  segment 
downrange.  From  the  differential  form  of  the  Biot-Savart  law 
(Ref.  9)  the  magnetic  induction  at  (r',  z'  ) attributable  to  the  n^1  seg- 
ment is 


B 

n 


Ih 


(r'  z') 


f lo 


1 x'f' 

An 


L 


Q6- 


(n2]3/2 


(19) 


where  T"is  the  magnitude  of  the  current  segments.  Allowing  for  re- 
versals in  current  from  segment  to  segment,  the  induction  produced 
by  all  transposed  segments  then  is 


®tr  p 


I 

til  X T OO  , Q6. 

is s (-i)n+1  ; — 

An  n=  1 (n-l)^[(zj 


dz ' 

z')2  + (r  ")2] 3/2 


(20) 


which  upon  integration  yields 


uj  x?'  ~ 
0™P<r'.  ^ ~ " 


zj  - (n  — Da 


Vt  z j - (n  - 1)^12  + (r')2 

(21)  ' 


In  the  limiting  case  of  an  infinite  line  current,  i.  e.  , n = 1,  z - */ 2, 
and  a-*  »,  the  summation  reduces  to  a quantity  of  2 so  that,  as  it 
should  per  Ref.  10, 


0I N F^r  ’ z ) 


T A' 

Po  1 X f 
2nr' 


(22) 
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The  ratio  of  transposed  magnitude  to  nontransposed  magnitude  then  is 


A plot  of  this  ratio  truncated  at  n = 9 has  been  obtained  through  use  of 
a digital  computer'"'  and  is  shown  for  several  values  of  */r"  in  Fig.  8. 
From  this  it  should  he  clear  that  if  *-/b  is  15  to  20,  as  it  will  be,  then 
a traveling  projectile  will  see  the  current  segments  as  essentially 
infinite  in  length  80  percent  of  the  time,  The  remainder  of  the  time 
is  spent  near  transpositions  where  the  radial  and  axial  fields  of  the 
angular  current  segments  become  effective.  Again,  the  shell  sees 
alternating  polarities  of  these  at  relatively  high  frequencies  and 
therefore  experiences  a system  of  diamagnetic  forces  that  on  the 
average  act  repulsively  toward  the  guideway  centerline. 

In  view  of  the  above  arguments,  a reasonable  approximation  is 
that  the  conductive  shell  experiences  the  same  forces  as  if  it  were 
traveling  through  the  null -flux  field  of  four  infinite  line  currents  but 
was  diamagnetic  to  a degree  corresponding  to  the  currents  (and  fields) 
being  segmented  and  transposed.  Fortuitously,  again,  a treatment  by 
Smythe  (Ref.  8,  p.  309)  relative  to  the  placement  of  a cylinder  of 
magnetic  material  with  permeability  q adjacent  and  parallel  to  a 
single  such  current  I is  applicable.  With  the  radius  of  the  cylinder 
being  a and  the  distance  from  its  axis  to  the  line  current  being 
r = r^,  as  shown  in  Fig.  9a,  the  significant  results  are: 

1.  The  magnetic  vector  potential  in  the  region  outside 
the  cylinder  due  to  its  presence  is  the  same  as  if 
the  cylinder  were  replaced  by  a parallel  "image" 
current  I'  located  between  I and  the  axis  at  a dis- 
tance a^/ri  from  the  latter,  plus  a current  of  -I' 
along  the  axis  as  shown  in  Fig.  9b.  The  magni- 
tude of  the  image  currents  is  given  by 

] - v-Jv- 

I = I 

l + (24) 


jijj 

See  Appendix  A for  a listing  of  Fortran  programs. 
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2.  The  vector  potential  inside  the  cylinder  is  the  same 
as  if  the  cylinder  were  removed  and  I were  replaced 
by  I",  where 


I " = — r-  I 

1 + f*  Jv  (25) 

The  latter  result  is  directly  comparable  to  the  shielding  effect 
of  a conductive,  cylindrical  shell  by  virtue  of  the  fact  that,  within  its 
interior,  the  vector  potential  of  an  exterior  line  current  I lying 
parallel  to  its  axis  is  the  same  as  if  the  shell  were  removed  and  I 
were  replaced  by  1"  where 


l"  = Si 

From  Eqs.  (25)  and  (26)  there  then  is  the  equivalence  that 


(26) 


1 + ^ (27) 

Upon  solving  this  for  p0/^  in  terms  of  S and  substituting  in  Eq.  (24) 
there  results  the  fact  that  when  there  is  effective  shielding  the  vector 
potential  outside  the  conductive  shell  due  to  parallel  line  currents  is 
the  same  as  if  the  shell  were  replaced  by  image  line  currents  I'  and 
-Fas  shown  in  Fig.  9b  where 


T = (S  - 1)  I (28) 

The  force  per  unit  length,  F,  on  the  cylinder  is  then  the  same  as 
the  sum  of  that  which  would  be  exerted  on  the  two  image  currents  by 
the  field  of  the  original  current,  i.  e.  , 


F = -r  X B,(0,  0)  + V X Bj(a2/rj,  0)  . (29 

From  Eq.  (22),  with  r'  being  appropriately  replaced,  the  magnetic 
induction  of  I at  the  points  of  interest  are 


Bj(0,  0)  = - 


2ttt ^ 


(30) 
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and 


13, {a"  'r,,  0) 


Ko1  rI 


(31) 


In  conjunction  with  appropriate  vector  identities  combination  of 
Eqs,  (28)  through  (30)  yields  finally,  for  the  force  per  unit  length  on 
the  cylinder 


_ fi  1 2d2($  - Dr* 

F = -- 

^r,(rf  - .,2)  (32) 

Because  S is  relatively  near  zero  for  the  diamagnetic  case,  it  is 
seen  that  the  force  is  directed  opposite  to  r,  i.e.  , away  from  the 
original  current  I.  A plot  of  F versus  ri  is  given  in  Fig.  10. 


For  the  situation  of  interest  herein,  i.  e.  , four  line  currents  of 
equal  magnitude  surrounding  the  cylinder,  the  image  currents  are  as 
shown  in  Fig.  11.  Note  that,  because  of  equal  magnitude  Ip  and 
equal  quantities  of  go  and  return  currents,  there  is  complete  cancel- 
lation of  images  located  on  the  cylinder  axis.  However,  for  each  line 
current  In  at  distance  rn  from  the  axis  there  is  a corresponding  image 
I'n'  that,  as  before,  lies  a distance  a2/rn  away  from  the  axis  toward 
In.  With  radial  and  angular  position  of  the  cylinder  axis  being  denoted 
by  51  and  (3,  respective!}',  cartesian  coordinates  xn  and  yn  of  the  line 
currents  have  been  taken  to  be 


X 

X 


X 


3 


-b 


X 


4 


0 


Yl  = 0 

y2  = b 
y3  = 0 


(33) 


By  triangulation,  coordinates  of  the  images  then  become 


^ a2(b  - 3?  co  n j3) 

*]  ' = 51  COS  & + „ .9  c 

1 5(2  + - 2J\b  cos 


<2  ' = S cos  f3  - 


Xg  - = 51  cos  0 - 


j251  cos  /3 


+ b*  — 2Xb  sin  j8 

a-(b  + 51  cos  (3) 

S2  + b2  + 251b  cos  fi 


y i ' = 51  sin  /3  - 


;i"51  sin  /3 


51-  + h2  — 251b  cos  f3 


-*sin/3  + a2(b  ~ Sln^)  <34) 

5i2  + b2  - 251b  sm  /S 


y3  ' = 51  sin  — 


a251  sin  f3 


51-  + b-  + 251b  cos  j3 
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K cos  /S  — 


i25t  co  s j3 


SI2  + b2  + 2&b  sin  jS 


*4 


<•  = 5?  sin  /8  — 


a2(b 


K sin  /l) 


X2  + b2  + 2fRb  sin  /9 


(34) 


The  magnetic  force  per  unit  length  experienced  by  the  cylinder  is 
the  same  as  that  which  would  be  experienced  by  the  four  image  currents 
in  the  combined  fields  of  the  original  line  currents  and  can  be  written  as 


F 


M 


1 (-l)n  Ip  x yn^ 

y i 


Now  from  Ref.  11 


Bi  (v-yn^  = Vx  Ai(x>  y>l 


(35) 


(36) 


and 


where 


A;(x,  y) 


fVp 

277 


4 


2(-l)°ln  , 
n=  1 D 


%=  \/<*  - Xn)2  + (y  - ya)2 


(37) 


(38) 


Also,  for  this  two-dimensional  case  with  i and  j being  unit  vectors  in 
the  x and  y direction,  respectively. 


_ — <?azA  <5az  ^ 

v x a(x,  y>  = : - -jt  i 


(39) 


so  that 


B.(x  / y -)  = — i(-Dn 
n’^n'  277  n=  L 


yn'-yn 


( X'-  x )2  + (y  y )* 
d a * n 3 n 


{v  - *n)2  + (v- 


(40) 

Through  combination  of  this  result  with  Eqs.  (28)  and  (35)  there  re- 
sults for  force  per  unit  length  on  the  cylinder 
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f„  = _ i)  i (_!)“'  i(_«» 

n =]  n=  ] 


M 


y-n'- 


( xn  “ x„)  + (y„' — y ) 

u n J n ■ n 


(x  ' - x J + ^Yn'  - y J2 


(41) 


This  can  be  resolved  into  radial  and  angular  components  by  taking  dot 
products  with  unit  vectors  so  directed,  i.  e. , 


(42) 


and 


(43) 


where 


and 


3!  = i cos  j3  + j sin  /3 


ft  = i sin  /3  — j cos  /3 


The  resulting  components  are 


(44) 


(45) 


^ - # « - » j«*  ;H  <4a, 

( n — *n  + (yn-  - >□’  J 


and 


1 m/9 


- (S  - 1,  i (-,)»•  X (.1).  K'~  - <v- 

t <v-  V* +<>■•- V*  ] 


(47) 


These  equations,  in  conjunction  with  Eqs.  (33)  and  (34),  have  been 
normalized  with  respect  to  b and  solved  at  varying  positions  using  a 
digital  computer.  v The  results  for  various  values  of  a/b  are  shown 
in  Fig.  12  for  radial  forces  and  Fig.  13  for  angular  forces. 


A listing  of  the  Fortran  program  is  included  in  Appendix  A. 
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It  should  be  noted  from  Fig.  12  that  for  a/b  < 0. 2 the  radial 
forces  for  all  8 can  be  approximated  by  a single  straight  line  if 
$/ b < 0.3.  It  is  also  to  be  noted  from  Fig.  13  that  for  this  condition 
the  angular  forces  remain  zero  for  all  practical  purposes.  There- 
fore, for  this  condition  the  positional  forces  are  purely  radial  and  it 
is  possible  to  define,  for  later  use  in  development  of  flight  equations, 
a linear  positional  force  coefficient  K]vi  such  that 


v _ FM$b 

M 3t  (48) 

where  h denotes  the  length  of  the  conductive  shell.  When  using  the 
slope  of  the  straight  line  approximation  on  the  normalized  curves 
this  becomes 


[£h(S  - l) 

K« ^ >•  ■l°P*  (49) 

In  summary,  the  slope  of  these  normalized  values  is  shown  in  Fig.  14 
as  a function  of  the  ratio  a/b.  Use  of  these  results  is,  of  course, 
strictly  valid  only  when  the  position  of  the  projectile  is  such  that  the 
ratio  JR/b  is  <0.3  and  the  ratio  a/b  is  <0.2.  With  discretion,  use 
beyond  these  limits  can  obviously  be  permitted. 

An  alternative  which  would  be  advantageous  in  this  regard  would 
be  to  gradually  spiral  the  field  modulation  currents  of  Fig.  1 in  lieu 
of  making  sudden  transpositions.  One  complete  spiral  each  4*  would 
give  the  needed  fundamental  frequency  for  effecting  diamagnetic  be- 
havior. Such  spiral  would  be  so  gradual  that  for  any  z the  currents 
and  their  field  would  again,  for  all  practical  purposes,  be  two- 
dimensional,  To  a traveling  projectile,  the  radial  forces  experienced 
would  tend  to  be  an  average  of  radial  forces  for  varying  from  0 to  2n. 
This  would  result  in  the  radial  forces  being  approximated  by  a single 
straight  line  for  larger  values  of  3l/b  than  would  be  the  case  for  the 
previous  transposed  system.  For  example,  in  the  case  depicted  in 
Fig.  12a  the  straight  line  approximation  for  the  spiraled  field  should 
be  valid  at  least  for  ft/b  = 0.  4. 

3.1.3  Damping  Force  Relationships 

Consider  now  the  interactions  which  occur  between  the  conduc- 
tive shell  and  the  field  of  the  nontransposed  currents  of  Fig.  2.  First 
of  all  it  is  to  be  noted  that  if  the  shell's  motion  is  strictly  parallel  to 
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the  z axis  it  will  sense  no  changes  in  the  magnetic  field  to  which  it  is 
exposed.  Hence,  once  equilibrium  penetration  is  established,  no  eddy 
currents  will  be  induced  and  there  will  be  no  interactions.  However, 
it  is  clear  that  lateral  motion  or  rotation  will  result  in  the  shell  being 
nonuniformly  exposed  to  changing  amounts  of  flux  so  that  eddy  currents 
will  be  induced.  These  are  to  be  utilized  for  the  purpose  of  damping 
oscillatory  flight  excursions  across  the  range  centerline.  In  this 
case,  as  stated  earlier,  tailored  flight  is  required  so  that  eddy  current 
power  will  be  essentially  real  {>98  percent  real,  see  Fig.  6).  Other- 
wise, shielding  will  occur  and  randomly  orientated  positional  forces 
will  be  introduced. 

As  the  first  step  toward  relating  damping  forces  to  other  param- 
eters, the  distribution  of  eddy  currents  in  a cylindrical  shell  arising 
out  of  relative  motion  with  respect  to  a single  parallel  line  current 
ID  will  be  found.  Superposition  will  then  be  used  to  obtain  the  solution 
for  the  case  of  four  such  line  currents  in  the  null -flux  arrangement  of 
interest.  From  this,  the  power  to  joulean  heating  of  the  shell  will  be 
obtained  and  finally  related  to  retardation  or  damping  of  its  radial 
and  angular  motion. 

Now  consider  the  position  of  a streamline  on  the  surface  of  the 
cylinder  to  be  specified  by  polar  coordinates  a,  and  the  position  of 
the  line  current  to  be  specified  by  ri,  6\,  both  as  shown  in  Fig.  15. 
Also,  consider  the  end  surface  of  the  cylinder  to  be  highly  conductive 
so  that  each  is  isopotential  at  $q(t)  and  $2(t),  respectively,  as  shown. 
Next,  insofar  as  surface  conditions  are  concerned,  it  is  to  be  recalled 
from  Fig.  6 that  for  power  > 9 8 -percent  real  the  magnetic  vector 
potential  of  the  eddy  currents  is  relatively  small  in  comparison  to  that 
of  the_inducing  field.  Moreover,  the  phase  angle  e between  the  inducing 
field  Aj,  and  the  field  A'  of  the  eddy  currents  is  such  that  their  vector 
sum  differs  little  in  magnitude  from  that  of  the  inducing  field  Aj  alone. 
Hence,  for  the  situation  of  interest,  only  the  latter  need  be  considered 
insofar  as  solving  for  eddy  current  magnitudes  and  distribution  is  con- 
cerned. 

Returning  to  Fig.  15,  the  total  electric  field  along  any  z-directed 
streamline  of  the  cylinder's  surface  is  that  induced  along  this  line  plus 
that  caused  by  charge  accumulations  and  can  be  written  for  the  two 
dimensional  case  as 


As  would  result  from  positioned  forces  alone. 


23 


AEDC-TR-76-149 


— 5A£.(a,  0,  t)  ^ 

E(a,  0,  t)  -g + [0  (0  _ 4>,(t)]r 

3t  2 J h (50) 

so  that  from  Ohm's  Law  the  current  density  is 


(51) 

where  ? is  area  resistivity,  h is  the  cylinder  length,  and  AD  is  the 
magnetic  vector  potential  of  Id.  From  Ref.  12,  the  vector  potential 
can  be  expressed  as 

An(a,  0,  t)  = — ^ In  yjx\  + a2  - 2ar,  cos  (0  - 0j)  + Cj 


i(a,  0,  t)  = ^j- 


<?AQ(a,0,t) 

dl 


*>  (01 1 


with  the  argument  of  In  being  the  distance  from  the  streamline  of 
interest  to  the  line  current  and  with  Ci  being  an  arbitrarily  large 
constant.  Upon  substitution  of  this  into  Eq.  (51)  and  differentiation 
with  respect  to  tipie  there  results 


i( a,  0,  t)  = 


t*ol  D (trl  - a coa  (0  — d j)l  dr  j/ dt  — arj  sin  (0  - (?j)d0j/dt| 


277? 


+ Tj  - 2 arj  cos  (0  ~ 6 j) 


<&2(t)  - 4>j(t) 
Ch 


(53) 

In  conjunction  with  the  fact  that  i has  only  z conponents,  an  obvious 
boundary  condition  is  that  there  can  be  no  net  flow  of  current  past 
the  cylinder  ends,  i.  e. , 


O 


(54) 


After  substituting  for  I (a,  0,t)  from  Eq.  (53),  integration  of  Eq.  (54) 
leads  to  the  result  that 


4«2(t)  - Ot(t)  fi0iD  drj 

h 2irrj  dt  j"  e 

Upon  returning  this  result  to  Eq.  (53)  and  manipulating,  the  distri- 
bution of  eddy  currents  attributable  to  one  line  current  becomes 
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i{a,  </>,  t)  - 


^£qFq  { arj  cos  (t£  — 0^)  - az  dr^ 

^n<'  a2  + r2  — 2 arj  cos  (<£  — 0j)]  dt 

ar | sin  (<f>  — 6 () 


(56) 


+ r2  — 2 aq  cos  ((/>  — $) 


dt 


In  generalized  form  for  the  h*-*1  of  several  line  currents  this  be- 
comes 


ia(a,  i)  = 


f^D 

2rr 


r cos  {<£  — 9)  — 

n - n 


f'Js2  + rn  “ 2arn  cos  ^ ~ ^ 


dt 


ar  sin  (<£  — (9  ) 

u II 


(57) 


+ — 2arn  c°3  ^ — 


dfl 

n 

IT 


For  the  situation  of  interest  which  is  four  line  currents  of  equal 
magnitudes  and  alternate  directions  as  shown  in  Fig.  16,  the  total 
current  density  can  then  be  written  as 


— i s o^d  ^ ( nn_ i 

*’£*’*'  t]  = 2^"  n=l  1.  r -2  . .2 


arR  cos  ((/>  — — 


ar  sin  (<A  — (9  ) 
n.  ' n 


rn[a‘*  + rn  - 2arn  co*  (<£  ” M 
d0_) 


d^n 

dt 


(58) 


+ r2  - 2 ar  cos  (<p  - $ ) 
n n “ n 


dt 


Again,  as  indicated,  the  Cartesian  coordinates  of  the  line  currents 
are  taken  to  be  <b,0),  (0,b),  (-b,0),  and  (0,  -b).  With  the  cylinder's 
radial  and  angular  positions  being  $ and  0,  there  then  exists  in 
addition  to  Eq.  (58)  the  generalized  geometric  relationships  that 


rQ  = \/4)2  + S2  ->  2‘Rb  cos  ■ - 7T  - /3) 


(59) 


and 


a n-l  . 

= rr  + sin 

n 2 


[r  sin  fr1  - 


(60) 


Differentiation  of  these  with  respect  to  time  leads  to  the  further 
generalized  relationships  that 


d/3 

dt 


dr 

n 

iR  — b CO  3 

(r’-'f 

dfR 

$b  sin 

1 

I 

k 

M 1 

ih 

"dt"  = 

rn 

dt 

_ 

r-  J 

(61) 
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and 


(62) 


Now  recall  that  the  radial  and  angular  components  of  the  shell's 
velocity  are 


and 


(63) 


'P 


= %T- 
At 


(64) 


Upon  substituting  these  into  Eqs.  (61)  and  (62),  the  results  into 

Eq.  (58),  and  regrouping,  the  total  eddy  current  distribution  resulting 

from  the  four  currents  becomes 


*<]*(  4> ’ 


*Vd 

2*C 


S(-l)‘^1[?n(a,b,  r,  % 


+ § (a,  b,  r,  % p,  <f>,  0) 


(65) 


where 

^ [arncos<0  — — a^l[5{  — bcos(-^-rr  - jS)] 

r2[r2  + a2  — 2ar  cos  (0  — d )] 
n n n “ n 

ar  sin  (<£  ~ B ) sia  (— — /3) 

it 

[r2  + a2  - 2arn  cos  (tf>  - 0n>][r2  — 5?2  sin2  17  ~ 

Sa  sin  (0  — sin  (— n— /3)[  $ — b cos  (^—  n — /9)] 

^[r2  + a2  — 2 arn  cos  — 0n)][r2  — JR2  sin2  — j8)]^ 


(66) 
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and 


[ar^  cos  (<p  — f?n)  — a2]  [— b ^in  ( n — /3)] 

'n^n  + a2  “ 2arn  cos({/>  " M 

[arn  sin  W ~ 6)n)]COfi(^'7  ~ P* 

+ [,2  + a2  - 2arn  cos  (0  - (MKr2  - X2  sin2  l^n-  /S)]’* 


§n  = 


Sab  sin  (<£  — 0 ) .sin ^ l 77  — /3 1 


rJL'l  + a2  - 2 arn  cos  (0  - 0n>][r2  - S2  sin2 


(67) 


Now  the  power  dissipated  as  joulean  heat  per  unit  cylinder  length 
is 


P = ac  / [TT(a,  d>.  t)]2d$  (68) 

O 

If  V^g  is  zero,  that  power  per  unit  length  attributable  to 
becomes 


* pKv  = ^ i(-r1J>b,r,iMiv^ 

4i72C  o n=l  (69) 

But  this  power  is  derived  from  mechanical  motion  where  the  incre- 
mental work  per  unit  length  is 


= -Fsvs 


(70) 


with  Fflfyl  being  the  implied  mechanical  force  acting  radially  per  unit 
length  against  a retarding  or  damping  force  , The  instantaneous 
power  or  rate  of  doing  mechanical  work  per  unit  length  is  then 


(71) 


Now,  for  later  convience,  radial  and  angular  magnetic  damping 
coefficients  are  defined  such  that 

kMdS  = r5th  ^ (72) 
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and 


K\ID0  “ '*£h  V1  (73) 

Equating  Eqs.  (69)  and  (71)  and  combining  with  Eq.  (72)  results 
finally  in  ■ 


KMd5(  = ~ 


9 O 

3)1  'h 


A /3,  c6.  0ij2  <10 


(74) 


By  setting  VR  equal  to  zero,  the  angular  damping  constant  is  obtained 
in  a similar  way  with  the  result  that 


K 


ah  2 


VI  D/3 


4>72  C 


° 1 U)-ISn(a,b.r,J!,/3.^0)]2^ 

O Ic=l  ^ 


(75) 


These  equations  have  been  normalized  with  respect  to  "b"  and 
solved  with  a digital  computer-  for  the  same  positions  as  before  in  the 
treatment  of  diamagnetic  forces.  Surprisingly,  for  any  given  position, 
KMDfR  and  KMD/3  are  equal.  Normalized  results  for  both  are  shown 
in  Fig.  17.  It  is  to  be  noted  that  within  the  previous  limits  established 
for  linearization  of  positioned  forces,  i.e.,  a/b  < 0, 2 and  ji/b  <0.3, 
the  damping  coefficients  are  essentially  constants  that  are  independent 
of  8 as  required  for  a linear  damping  term  in  the  equations  of  motion 
to  be  developed  later.  In  either  the  positional  or  damping  situations 
a choice  of  larger  limits  would  not  be  unrealistic.  In  summary,  their 
normalized  value  is  plotted  versus  the  ratio  a/b  in  Fig,  18. 


3.2  FLIGHT  DYNAMICS 

For  simplicity,  the  following  development  will  assume  flight  to  be 
confined  to  a v-z  plane  of  symmetry  which  passes  through  both  the 
projectile  axis  and  the  guideway  centerline  (z  axis).  This  assumption 
is  justified  on  the  basis  that  the  projectile  is  to  be  a body  of  rotation 
and  the  guideway  design  is  to  be  such  that  magnetic  positional  forces 
are  almost  totally  radial.  Any  angular  components  are  sufficiently 
small  that  their  effect  on  dispersion  of  the  projectile  from  the  guide- 
way centerline  should  be  negligible.  A means  of  visualizing  their 
effect  is  to  imagine  relatively  slow  rotation  of  the  plane  of  symmetry 
about  the  z axis. 

^A  listing  of  the  Fortran  program  is  included  in  Appendix  A. 
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Within  the  plane  of  symmetry  y will  boused  to  denote  the  angle 
between  the  z axis  and  the  velocity  vector  V of  the  projectile's  Cg 
{center  of  gravity),  a will  denote  the  projectile’s  angle  of  attack  with 
respect  to  V,  and  v will  denote  the  dispersion  of  the  Cg  from  the 
z axis  {all  shown  in  Fig.  19).  Still  within  this  plane,  the  aerodynamic 
lift  and  drag  forces  will  be  considered  to  act  through  an  assumed 
center  of  aerodynamic  pressure  cp  a distance  i aft  of  cg,  and  the 
magnetic  positional  and  damping  forces  will  be  considered  to  act 
through  an  assumed  center  of  magnetic  pressure  c^[  a distance  fj/j 
forward  of  Cg, 

Recall  also  from  Figs.  12,  13,  and  17  that  the  design  can  be 
such  that  the  magnetic  positional  forces  will  be  a linear  function  of 
distance  from  the  z axis  with  the  magnetic  damping  forces  being  a 
linear  function  of  the  rate  this  distance  changes.  Additionally,  the 
allowable  angles  of  attack  are  to  be  small  so  that  aerodynamic  lift 
and  damping  can  be  assumed  linear  functions  thereof  with  drag  an 
independent  linear  function  of  velocity  alone.  Then  at  fixed  velocity 
within  the  anticipated  operating  conditions  the  equations  of  motion 
can  be  assumed  to  be  linear  and  will  be  so  treated. 

3.2.1  Equations  of  Motion 

Referring  again  to  Fig.  19,  the  sum  of  y forces  acting  on  the 
projectile  must  be  zero  so  that  for  small  a and  y one  can  write 


M — ^ - qSbCLaa  ~ + ^M(a  + + qSbCDO>'  " KMD^d9  dt>+  - 0 

dt2 


where  previously  unmentioned  parameters  M and  are  projectile 
mass  and  base  area  respectively,  q is  dynamic  pressure  (IpV^  with 
p being  density  of  the  working  gas),  and  Cl  a an<3  Cjyo  are  t*ie  aer0“ 
dynamic  lift  and  drag  coefficients,  respectively. 


Additionally,  the  sum  of  moments  about  cg  must  be  zero  so  that 
one  can  also  write 

♦ 9V<L  + CD0)“  - 2iVbcw  r'l 

dt2 

- + tM(a  +■  y)l  - kMD£MP<Vtlt)  + £M((k/dt)]  = 0 (77) 
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where  J and  rb  are  projectile  moment  of  inertia  and  base  radius, 
respectively,  and  is  the  aerodynamic  damping  coefficient. 


The  angle  y and  its  derivative  can  be  eliminated  from  Eqs.  {76) 
and  (77)  by  the  approximation  for  small  angles  that 


Y = 


(78) 


Upon  doing  this,  the  Laplace  transformation  of  the  resulting  equations 
has  made  it  possible  to  combine  and  solve  algebraically  for  both  y and 
a in  terms  of  the  complex  quantity  s with  the  result  that 

KD6  + VOyJ  *3  + [(DSV  - D2  + D0)yo  + (D6  + V)y0]s2 

- KDA  - °4  + DlD6>*o  - D6V«c  - DS  Vflo  + D6»°  - 

+ [{D4D5  - DjDg)^o  - <D4D6  + D pB)ao  + (D4*o  + Dgyo  + D8Vao)] 

«(s)  = 

l(D6  + V)*4  + [D0  + D5v  - DJs3  + [(D4  - D2D5  - DyV  - DjDgls2 
+ [D4D5  + D2D7  - DA  - d,d8]9  - [D3D8  + D4n7] 

(79) 

and 

[(Dg  4-  V)a0]ss  + [y0  - D^zo  + Vdo  + D54o  + D$Vao  - D?yo]  s2 
+ [(Djyo  _ DjD6ao  + Dgy0  + D^o  - + D5Vao  - D?  y0-D7Vao]s 

a(s)  = ~ [D3D6ao  " D1D7Vq  “ P3*o  ~ D3*Vo  + P7^o  ~ D2D7ao  + 

t(D6  + V)3s4  + [Dg  + Dsv  - DjJs3  + [(D4  - D2D5  - D7V  - DjD6]s2  “ 

+ + D2D-  - D3D6  - DjDgls  ~ [D3Dg  + D4D?] 

(80) 

where  the  subscript  "o"  denotes  the  initial  condition,  the  superscripts 
* and  • • denote,  respectively,  the  first  and  second  derivatives  with 
respect  to  time  and 


Di  - + kmd?mv^  1 

D5  = [qSbCD0  - - K^VlM-W-l 

D2  = [2^brbCnfi  + KMDgMV^J  1 

°6  = 

(81) 

d3  = v^vj-1 

d7  = kmm-! 

D4  = VCqSgCfC^  + CDq)-  KmBm1  j 

Dg  = [qSbCu  + 
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Note  that  the  denominators  of  Eqs.  (79)  and  (80)  are  identical  quartics 
and  that  the  numerators  are  cubics  whose  difference  lies  in  the  coeffi- 
cients of  the  powers  of  s.  Letting  the  roots  of  the  quartic  be  Hj,  R2, 
R3,  and  R4,  a general  expression  for  Eqs.  (79)  or  (80)  is 

,,  . as^  + + cs  + <1 

f(.s)  = 

9(S-R1){S-R2HS-R3)(S-R4)  (82) 


where,  for  brevity,  a,  c,  d,  and  » are  used  in  place  of  the  coeffi- 
cients. Upon 'application  oi  Heaviside's  Expansion  Theorem  (Ref.  13) 
an  equivalent  expression  in  terms  of  partial  fractions  is 

K.  K„  K,  K, 

«•>  - — * — + — * — 

s — R ^ s-Rg  s—  Rj  a— too  7 

where 

Kt,  = |(s-Rn)f(8)j^Rn  (84) 


i.  e.  , 


(s-  Rn)ia  Rp  + R„  + + dl 

9.(3  ~ Rj)(s  - r2)(s-  r 3)(s  - n4) 


(85) 


An  inverse  transformation  of  Eq.  (83)  yields  finally  that  in  the 
time  domain 


f(0  = 


S ,K„exP  (M 
0=1 


(86) 


Recall  that  this  is  a general  solution  which  is  applicable  to  solving 
Eqs.  (79)  or  (80)  for  tf(t)  or  cr(t). 

Unfortunately,  the  solution  does  not  yield  a general  expression 
for  the  roots  of  the  quartic.  It  has  been  necessary  therefore  to  re- 
sort to  numerical  techniques  and  the  digital  computer  as  a means  of 
obtaining  their  value  for  each  set  of  chosen  parameters.  At  the  same 
time  the  computer  has  also  been  used  to  determine  the  values  of  Kj, 
K2,  K3,  K4,  and  finally  to  yield  a timewise  solution  for  both  v and  a. 
A listing  of  the  Fortran  program  is  included  in  Appendix  A, 
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3.2.2  Discussion  of  Roots 

In  general,  as  might  be  expected  from  the  physics  of  the  problem, 
the  four  roots  have  been  found  to  consist  generally  of  two  pairs  of 
complex  conjugates.  One  pair  represents  an  oscillation  that  exists 
primarily  because  of  aerodynamic  forces.  The  second  pair,  whose 
value  is  also  influenced  greatly  by  aerodynamics,  arises  when  mag- 
netic guidance  is  added  which,  to  a certain  extent,  modifies  the  first 
pair.  The  interrelation  is  involved  and  a clear  separation  has  not 
been  made. 

A design  methodology  which 'would  allow  direct  determination  of 
unknown  parameters  from  initially  specified  roots  might  possibly  be 
developed  in  the  future.  The  means  of  achieving  this  would  be  through 
use  of  the  roots  to  form  a quartic  as  in  the  denominator  of  Eqs.  (79) 
and  (80),  The  methodology  would  consist  of  a systematic  means  of 
adjusting  parameters  to  yield  equivalent  coefficients  for  all  powers 
of  s. 


4.0  RANGE  "G"  CORRELATION  AND  IMPLEMENTATION 

The  following  preliminary  adaptation  of  the  concept  to  the  VKF 
Range  G is  for  demonstration  purposes  only,  is  limited  in  scope,  and 
is  not  to  be  taken  as  final. 

4.1  BASIC  PHYSICAL  PARAMETERS 

To  minimize  scaling  to  actual  reentry  situations,  it  is  generally 
preferred  that  tests  in  aeroballistics  ranges  be  performed  on  the 
largest  possible  scale.  Toward  this  end  a base  diameter  of 
5.08  x 10“2  m will  be  assumed  for  the  projectile  proposed  herein. 
This  is  approximately  the  largest  that  can  be  accomodated  by  the 
6. 35  x 10-2  m bore  of  the  Range  G launcher  in  view  of  the  fact  that 
during  launch  the  projectile  is  to  be  surrounded  by  a sabot  which 
subsequently  separates.  The  cylindrical  portion  of  the  projectile, 
and  hence  the  diameter  of  the  conductive  shell,  will  be  assumed  to 
be  one-half  the  base  diameter,  i.  e.  , 2.  54  x 10-2  m.  The  length  of 
the  shell  will  then  be  taken  as  four  times  its  diameter,  i.  e.  , 

1.016  x lO--*-  m.  This  is  a reasonable  compromise  between  a desire 
on  the  one  hand  for  a small  length-to-diameter  ratio  to  withstand 
launch  loads,  and  a desire  on  the  other  hand  for  a large  length-to- 
diameter  ratio  so  as  to  be  compatible  with  the  previously  developed 
theory. 
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The  launch  package  as  presently  envisioned  from  the  above  is 
depicted  in  Fig.  20.  As  indicated,  the  magnetic  forces  are  assumed 
to  act  at  the  midpoint  of  the  copper  shell.  " Also,  as  indicated,  it  is 
presumed  that  through  use  of  appropriate  flare  angle,  flare  length, 
deployment  of  material,  etc.  , the  projectile's  center  of  gravity  will 
be  placed  a suitable  distance  H m behind  and  the  center  of  aerodynamic 
pressure  a suitable  distance  1 still  further  behind  the  shell's  mid- 
point. For  the  purposes  of  this  study,  these  distances  or  moment 
arms  will  be  taken  to  be 


and 


£ = 9.525  x 10  3 m 


(87) 


eM  = 1.905  x I0'2m 

Perhaps  surprisingly,  the  influencing  factor  insofar  as  choice  of 
guideway  parameters  is  concerned  derives  from  the  sabot.  This  is 
because  the  first  30  m of  postlaunch  flight  is  through  a so  called 
"blast  chamber"  within  which  sabot  separation  and  destruction  occurs. 
Because  it  would  otherwise  interfere  with  this  event,  as  well  as  sus- 
tain damage,  the  logical  beginning  of  the  guideway  should  be  immedi- 
ately downstream  of  the  blast  chamber  exit.  At  this  point,  based  on 
almost  ten  years  of  G-range  history,  a 5.08  x 10'2  m dispersion  of 
the  projectile  away  from  the  guideway  centerline  must  be  assumed. 
Upon  adding  to  this  the  base  radius  of  the  projectile  and  an  allowance 
for  the  conductors  needed  to  carry  guideway  currents,  the  minimum 
allowable  choice  for  the  guideway  dimension  b {see  Fig.  11)  is 
8.  89  x 10-2  m.  This  value  will  be  utilized  herein.  Next,  from  Fig.  8 
it  is  seen  that  a good  choice  for  the  length  4 of  guideway  transpositions 
is  17  times  this  value  of  b or,  rounded  off,  1.5  m. 

Returning  now  to  the  projectile  itself,  the  design  velocity  is  taken 
to  be  6,000  m/sec.  The  effective  fundamental  angular  velocity  of  the 
spacewise  magnetic  field  modulations  then  becomes 

co  = 2ffV(2*)  ' = 12,560  radian  s ' sec  (88) 


This,  in  conjunction  with  the  correlation  parameter  of  Fig.  6,  enables 
determination  of  required  shell  area  resistivity. 
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As  seen  in  Fig.  6,  a reasonable  compromise  insofar  as  concerns 
diamagnetic  behavior  of  the  conductive  shell  is  where  the  eddy  currents 
are  98-percent  reactive  and  where  S has  a value  of  0.22.  At  this 
point  the  correlation  parameter  has  a value  of  4.  5,  i.e. , 

H0  a at  (2  n<? ) 1 = 4.5  (89) 

Now,  from  Ref.  8 (p.  299),  the  vector  potential  at  polar  coordinates 

r,  9 attributable  to  a line  current  lying  parallel  to  the  z axis  at 

ro,  60  can  be  expressed  in  terms  of  circular  harmonics  as  (see  Fig.  15) 


AZM)  = ^[ln  rc 

7 '2n  c 


— X ~ (r,r0)n  (cos  n 0o  cos  nd  + sin  n£?0  sin  n 6] 
n=  1 


(90) 


Also,  as  stated  earlier,  because  the  shell  is  of  finite  length,  the  net 
flow  of  current  past  its  ends  must  be  zero.  Hence,  for  the  postulated 
two-dimensional  situation,  any  induced  electric  field  attributable  to  the 
time  derivative  of  the  term  "in  r0M  of  Eq.  (90)  will  be  ineffective  in 
producing  eddy  currents.  All  which  exist  must  therefore  be  attribu- 
table to  the  time  derivative  of  the  indicated  summation.  Upon  dif- 
ferentiation of  the  latter  it  becomes  clear  from  inspection  that  if 
r0  is  somewhat  larger  than  r the  n = 1 term  dominates.  Hence,  by 
virtue  of  the  fact  that  during  flight  the  distance  between  the  projectile 
axis  and  any  of  the  line  current  segments  of  Figs.  1 or  3 is  to  be 
maintained  relatively  large  ,£  as  compared  to  all  r < a,  it  is  seen  that 
the  correlation  parameter  need  only  be  considered  for  the  value  n = 1. 
Upon  substituting  this  value  for  n along  with  the  previously  established 
values  for  a and  u into  Eq.  (89)  and  solving  for  area  resistivity  there 
results 

^ = 2.227  x 10 '5  ohms  (91) 


Assuming  the  shell  to  be  oxygen-free,  high  conductivity  copper  with  a 
volume  resistivity  of  p = 1.  724  x 10”^  ohm-m,  the  implied  shell 
thickness  is 


w =-£  = 7.74  x Hr4m 


(92) 


For  previously  stated  reasons,  the  thickness  of  the  end  walls  will  be 
taken  at  twice  this  value. 


All  dimensions  are  now  defined  that  are  needed  for  calculating 
the  mass  and  moment  of  inertia  of  that  portion  of  the  shell  (with  its 


*From  three  to  seven  times  greater. 
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Lexan®  interior)  which  lies  forward  of  the  projectile's  Cg.  The  total 
mass  and  moment  of  inertia  of  the  projectile  are  simply  twice  the 
result  plus  an  allowance  for  a nose  piece.  Upon  using  8.  89  x 10^  and 
1. 19  x 103  kg/m3  for  the  density  of  copper  and  Lexan,  respectively, 
there  results  for  the  projectile  totals 


M = 1.65  x 10  1 kg 


and 


J = 3.0  x 10‘4  kg  - m2 


03) 


4 2 GUIDEWAY  CURRENTS 

For  the  purpose  of  this  preliminary  correlation  to  the  G range 
the  following  aerodynamic  coefficients  are  believed  to  be  adequate 
and  will  be  utilized  throughout 


C 


DO 


0.3 


1.1 


-2.0 


(94) 


Only  flights  in  air  at  atmospheric  pressure  and  at  one -half  atmos- 
pheric pressure  will  be  considered.  For  these,  the  corresponding 
densities  are  taken  to  be  1. 1767  and  0.  5883  kg/m3,  respectively. 


4.2.1  Current  Magnitude 

For  the  particular  dimensions,  shielding  ratio,  and  area 
resistivity  established  above,  is  is  found  through  use  of  Figs.  14 
and  18  that 

Km  = -8.0  x 10-7I2 

and 


kmd 


-3.3  x 10“ 10 1 q 


(95) 
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Upon  insertion  of  proper  current  levels  the  resulting  values  of  Km 
and  Kjyjj-),  in  conjunction  with  values  of  the  other  parameters  as 
have  now  been  established,  form  the  input  data  for  the  computer 
programs  mentioned  in  Sect.  3.  2. 1.  Dispersion  (*  ) of  the  model 
about  the  range  centerline  and  its  angle  of  attack  (a)  as  functions 
of  downstream  position  as  well  as  roots  of  the  characteristic 
equation  of  Eqs.  (79)  and  (80)  can  be  computed. 

Being  primarily  influenced  by  the  current  rating  of  the  Tunnel  "F" 
generators  (Ref.  7)  the  authors  have  chosen  105  amp  as  being  an 
appropriate  Ip  and  Ip  magnitude  for  demonstration  of  the  concept. 

For  such  magnitudes  Eq.  (95)  yields 


and 


K\,  = -8  x 103N/m  (96) 

KMd  - -3.3  N - sec/m 


Initial  conditions  for  the  computed  examples  to  be  presented,  i.  e. , *Q, 
Q?OJand  v0,  are  in  all  cases  taken  to  be  5.08  x lCT^m  and  1 x 10~2  and 
1 x 10'3  radians,  respectively.  Results  for  an  atmospheric  case  in  the 
existing  range  length  of  300  m are  depicted  in  Fig.  21  while  Fig.  22 
demonstrates  the  effect  of  reducing  density  to  that  for  one-half  of  atmos- 
pheric. The  behavior  of  y and  a as  depicted  in  the  figures  is  in  agree- 
ment with  the  desired  effect,  and  the  roots  are  in  accord  with  expec- 
tations based  on  the  theory  of  controls. 

In  Figs.  23  and  24  the  velocity  is  reduced  to  4,  500  m/sec  with  an 
attendant  shift  in  shielding  ratio  and  accordingly  a change  in  magnetic 
postional  force  constant  to 


Km  = 6.154  N/m  (97) 

kMD  not  significantly  altered  since  its  magniture  is  not  dependent 
on  velocity  through  the  transposed  (or  modulated)  field.  All  other 
variables  remain  as  in  the  previous  two  examples,  and  Fig.  23  depicts 
the  atmospheric  and  Fig.  24  the  one-half  atmospheric  case.  Again, 
results  are  gratifying  and  guidance  is  demonstrated. 

The  effect  of  varying  the  levels  of  ID  and  Ip  can  be  seen  in  Figs.  25 
and  26  where  an  atmospheric  and  a one -half  atmospheric  case  similar 
to  those  of  Figs.  21  and  22  have  been  computed  except  using  an 
increased  Ip  of  1.25  x 10  5 amp.  It  is  apparent  that  the  reduction  in 
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error  from  the  initial  value  for  both  y and  a is  greatest  with  this  modi- 
fication. Reasons  for  presenting  the  other  examples  in  which  Ip  and 
Ip  are  the  same  will  become  clear  in  a following  section. 

Figures  21  through  26  exhibit  significant  reduction  of  initial  error 
with  any  of  the  above  mentioned  sets  of  parameters,  i.  e. , use  of  these 
sets  of  parameters  results  in  guided  flight.  It  should  be  emphasized 
here  that  these  results  are  in  marked  conflict  with  any  unguided  case 
where  for  the  assumed  initial  conditions  y would  have  attained  a value 
of  0.  3 5m  at  a distance  of  300  m downrange. 

It  should  be  noted  here  that  the  distance  between  excursions  across 
the  range  centerline  is  a minimum  of  45  m which  compares  to  a 1.5-m 
transposition  length.  Therefore,  as  desired,  eddy  currents  induced 
by  the  nontransposed  field  have  a frequency  more  than  20  times  lower 
than  those  resulting  from  velocity  through  the  transposed  field  and 
are  > 98-percent  real,  as  desired. 

A final  point  that  should  be  made  is  that  the  range  length  is  also 
an  important  parameter.  Note  that  if  this  were  increased,  with  all 
other  parameters  remaining  fixed,  the  initial  error  would  be  further 
reduced.  To  illustrate,  if  the  error  is  reduced  to  one -fifth  its  original 
value  in  a 300-m  range  it  would  be  reduced  to  a twenty-fifth  of  its 
original  value  in  a 600-m  range.  Hence,  application  to  longer  ranges 
is  encouraging  from  this  viewpoint. 


4.2.2  Implementation 

Although  not  specifically  stated,  the  establishment  in  Section  4. 1 
of  the  guideway  dimension  b (see  Fig.  2)  as  8.87  x 10"^  m includes  an 
allowance  for  conductors  of  diameter  2,  54  x 10’^  m to  carry  the 
guideway  currents.  This  diameter,  in  conjunction  with  the  b dimension 
and  a guideway  length  of  300  m are  the  basis  of  the  calculated  induct- 
ance and  resistance  values  that  are  used  in  the  following  development. 


Because  of  the  transposition  all  tendencies  for  mutal  inductance 
between  the  Ip  and  Ip  currents  of  Fig.  1 cancel  when  overall  guide- 
way length  is  considered.  Hence,  with  reference  to  Figs.  1,  2,  and 
16,  the  energy  stored  inductively  in  the  basic  guideway  field  of  the 
Ip  currents  alone  can  be  written  as 


HK(ind) 


4 [Ha- 


lf j 


- 2L 


Lmnt  1-3^  hi 


(98) 
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where  Lself  l is  the  self  inductance  of  each  of  the  conductors  and 
•^mut  1-2  anc*  is  the  mutual  inductance  between  diagonal 

and  adjacent  conductors,  respectively.  Based  on  working  formulae 
3 and  7 of  Grover  (Ref.  14),  approximate  values  of  self  and  mutual 
inductance  for  the  aforementioned  dimensions  are 

Lse]f  i = 6-01  x 10-4  h 

Lmut  i_2  - 4-48  x 30  4h 

Lmui  1—3  = 4.27  x 10*4  h (99) 

Upon  substitution  into  Eq,  (98),  there  results  for  inductive  energy  in 
the  Ip  field 


fflBF(md)  = 2.64  X 10-4  I*  {100) 

Although  consideration  of  transpositions  would  result  in  a slightly 
higher  value,  the  effective  self  and  mutual  inductances  of  the  conductors 
for  the  Ip  currents  will  be  taken  to  be  equal  to  those  for  Ip.  Using 
nomenclature  of  Fig,  11,  then  the  inductive  energy  of  the  transposed 
field  is 


®TF(ind>  - 2 64  x 10'4  Ip  (101) 

Total  inductive  guideway  energy  is  then 

®GW(ind)  = ®BF(ind)  + ®TF(ind>  = 2 64  x 10"*(!D  + (102) 

r 

For  Ip  and  Ip  of  10  amp  as  used  in  4.2.1  above  the  inductive  energy 
in  the  guideway  field  would  be 

^GW(ind)  = 5 28  x iq6J  (103) 

This  field  energy  level  coupled  with  current  levels  of  105  amp  and  the 
fact  that  the  guideway  need  only  be  energized  for  approximately  0. 1 sec* 
suggests  the  utilization  of  the  nearby  flywheel/ inductor  energy  store 
power  supply  for  the  VKF  Tunnel  F (Ref.  7).  Approximately  4 x 108  J 
is  available  from  its  flywheel- driven,  acyclic  generators  at  10 6 amp 
and  90  v,  or  at  5 x 10^  amp  and  180  v;  or  from  its  storage  inductor. 


',cThe  time  required  for  a projectile  to  travel  the  guideway  length, 
plus  contingency  time. 
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108  J is  available  at  106  amp  up  to  2 x 10^  v.  However,  unless  section- 
alization  is  utilized  as  will  be  discussed  below,  or  unless  conductors 
are  utilized  at  cryogenic  temperatures,  and  space  requirements  for 
thermal  insulation  discourages  this,  overall  guideway  resistance  will 
not  allow  either  approach.  If  it  is  of  copper  at  room  temperature, 
resistance  of  a single  conductor  of  2.  54  x 10 m diameter  and  300-m 
length  is 

n (300  m)(l.7l  x 10  ® ohm  — m) 

(rr.‘  4)  (2.54  X 10'2  m>2 

- 1.012  x 10" 2 ohm  (104) 


The  current  that  can  be  driven  through  two  such  conductors  in  series 
(one  go  and  one  return)  at  maximum  generator  voltage  of  130  v is  thus 


In  = 


180 

2(1.012  x 10'2) 


= 8,893  amp 


(105) 


This  is  inadequate  by  an  order  of  magnitude.  Moreover,  if  one  resorts 
to  utilization  of  the  storage  inductor  (assuming  its  current  could  be 
commutated  through  the  guideway  with  reasonable  efficiency)  the  energy 
to  joulean  heating  of  the  eight  guideway  conductors  by  a 10 b amp, 

10-1  sec  pulse  through  each  would  be 


®GW(rraD  = + IpKResHAt) 

= 4[(105)2  + (105)21(1.012  x lO^KlO-1) 

= 8 x 107  J (106) 

This  is  approximately  the  full  rating  of  the  storage  inductor. 
Considering  that  vast  amounts  of  energy  must  inescapably  be  lost  in 
commutating  current  through  the  guideway  and  that  the  delivered 
current  pulse  would  be  a rapidly  decaying  exponential,  this  approach 
is  unsatisfactory.  It  is  noteworthy,  however,  that  calculations  show 
the  desired  pulse  would  lead  to  a guideway  temperature  rise  of  only 
20°C. 


One  viable  alternative  would  be  to  procure  a large  flywheel  - 
driven  generator  whose  energy,  voltage,  and  current  ratings  would 
be  matched  to  the  guideway.  From  Eqs.  (101)  and  (106)  it  appears 
that  a total  flywheel  energy  on  the  order  of  the  6,  84  x 10 8 joule  total 
for  Tunnel  F would  be  adequate.  As  mentioned  above,  an  attractive 
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possibility  for  Range  G which  would  allow  utilization  of  the  Tunnel  F 
power  supply  is  offered  through  guideway  sectionalization.  Each 
section  would  be  powered  independently  of  all  others,  and  then  only 
for  the  period  needed  for  the  projectile  to  pass.  This  might  be 
achieved  by  explosively  or  perhaps  pneumatically  actuated  switches 
which  would  appropriately  connect  and  disconnect  each  section  to  two 
large,  range-length  conductors*  that  would  be  energized  by  the 
Tunnel  F generators.  Minimal  switch  developmental  effort  should 
be  required  because  of  the  low  operating  voltage  (180  v)  and  because 
the  guideway  inductance  to  resistance  ratio  is  relatively  low. 

As  still  another  alternative,  relatively  small  flywheel -driven 
generators  could  be  stationed  at  intervals  along  the  wayside  for 
powering  each  guideway  section  at  the  appropriate  time.  Or,  alter- 
natively, banks  of  energy  storage  capacitors  could  be  utilized.  It 
is  to  be  expected,  however,  that  the  rotating  machinery  approach 
would  be  the  most  economical,  more  trouble  free,  and  with  an 
expected  life  that  would  be  significantly  longer. 


Returning  again  to  Figs.  1 and  3,  attention  is  invited  to  the  fact  that 
if  Ip  and  Ip  are  equal  in  magnitude,  the  basic  guideway  field  is  canceled 
by  alternate  transpositions  and  is  doubled  in  value  by  the  ones  between. 
This  suggests  the  possibility  of  utilizing  currents  (and  their  conductors 
or  coils)  only  at  alternate  intervals  as  shown  in  Fig.  27.  To  make 
this  point  clear  is  the  reason  for  computing  flight  paths  in  Sec.  4,  2.  1 
for  equal  Ip  and  Ip  despite  the  fact  that  improved  results  were 
achieved  with  increased  Ip.  This  means  an  increased  level  for  the 
basic,  or  average  (nonmodulated  field)  is  desired.  Hopefully,  it 
can  be  seen  that  essentially  this  same  biasing  effect  can  be  achieved 
by  making  the  z spacing  between  the  coils  of  Fig.  27  less  than  the 
z dimension  of  the  coils  themselves.  An  obvious  bonus  feature  of 
such  spaced  coils  would  be  that,  while  in  the  intervals  between  coils, 
the  projectile  would  be  free  of  shock  interactions  and  there  would  be 
no  obstructions  insofar  as  wayside  photography  is  concerned 


"One  go  and  one  return  conductor,  each  having  an  area  on  the 
order  of  0.  1 m . 
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5.0  CONCLUDING  REMARKS 


An  electromagnetic -guideway/ guided-projectile  system  concept 
that  is  believed  suited  to  the  VKF  Range  G,  and  perhaps  numerous 
other  aeroballistic  ranges,  has  been  presented  and  developed  from 
a theoretical  viewpoint.  Basis  of  the  system  is  guided,  aerodynamic 
flight  via  motion-induced,  electrodynamic  interactions  between  way- 
side  magnetic  fields  that  form  a guideway,  and  conductive  material 
that  is  deployed  in  the  projectile  structure.  The  geometry  of  the 
fields  is  such  that  induction  decays  to  zero  with  proximity  to  the 
guideway  axis  so  that  interference  with  aerodynamic  flow  or  abla- 
tion near  a projectile's  stagnation  point  should  be  minimal.  Aside 
from  being  devoid  of  sliding  surface  contact  as  is  used  in  mechanical 
guideways,  a possible  attractive  feature  of  the  system  is  periods  of 
essentially  free  flight,  -with  no  obstructions  to  interfere  with  wayside 
photography. 

Quite  naturally,  system  theory,  and  hence  the  concept  itself, 
has  been  developed  on  the  most  elementary  basis  possible.  For 
example,  the  concept  as  presented  utilizes  conductive  material  in 
the  form  of  a cylindrical  shell  that  is  placed  principally  about  the 
projectile's  forebody.  It  should  be  possible,  and  from  an  aerody- 
namic viewpoint  perhaps  desirable,  to  achieve  the  necessary  elec- 
trodynamic interaction  with  the  guideway  if  a shell  of  conical  or  other 
shape  were  used.  However,  the  electrodynamic  theory  would  be  most 
difficult  to  develop  as  the  problem  would  no  longer  be  two-dimensional. 
Fortunately,  writh  the  cylinder  it  has  been  possible  to  develop  most  of 
that  needed  by  two  independent  means,  the  briefest  of  which  has  been 
presented.  Certain  idealizations  and  simplifying  assumptions  have, 
however,  been  made  during  the  course  of  the  development  which 
apply  to  both.  For  this  reason,  experiments  are  underway  that 
are  intended  to  confirm  theory  insofar  as  concerns  electrodynamic 
interactions.  A report  on  the  results  is  planned  for  the  near  future. 

No  experimentation  as  regards  aerodynamics  is  planned  as  adequate 
wind  tunnel  data  should  be  available  for  this  purpose. 

For  the  particular  case  of  Range  G the  possibility  of  energizing 
a guideway  from  the  nearby  Tunnel  F generators  is  optimistically 
offered,  provided  that  the  guideway  is  sectionalized  and  that  each 
section  be  powered  only  during  the  time  the  projectile  is  passing 
through,  plus  contingency.  Possible  alternate  methods  of  guideway 
energization  for  this  and  other  ranges  have  also  been  offered.  How- 
ever, in  retrospect,  it  is  possible  that  stronger  fields  with  corre- 
spondingly stiffer  flight  control  may  be  desired  in  lieu  of  these 


41 


AEDC-TR -76-149 


values  used  in  the  Range  G examples  of  Section  4.  2. 1.  One  reason  is 
that  ablating  and/or  eroding  nose  test  specimens  may  initially  have 
(or  tend  to  develop)  asymmetry  and  cause  the  projectile  to  become 
stubborn.  Another  is  that  greater  downrange  accuracy  may  be  needed 
for  capture.  This  could,  of  course,  be  improved  by  a smaller  initial 
error.  Such  reduction  additionally  would  allow  a smaller  guideway 
dimension  b,  and  hence,  from  Figs.  14  and  18,  greater  magnetic 
force  coefficients  for  given  currents.  Also,  with  sectionalization, 
there  is  the  possibility  of  utilizing  greater  currents  from  the  Tunnel 
F power  supply  which  would  lead  to  still  greater  force  coefficients. 

Still  another  reason  for  stiffer  control  is  to  achieve  a reduction 
in  the  maximum  angle  of  attack  assumed  by  the  projectile.  Quite 
obviously  this  would  be  achievable  with  lowered  initial  error.  But 
a reduction  could  also  be  achieved  by  virtue  of  the  fact  that  greater 
magnetic  forces  would  permit  lesser  reliance  on  the  development  of 
aerodynamic  forces.  In  turn,  lower  magnetic  moment  and  moment 
arms  would  be  required.  To  visualize  the  trend,  consider  the  case 
of  zero  where  the  magnetic  forces  would  attempt  to  correct  the 
projectiles  course,  but  with  zero  change  in  angle  of  attack.  Finally, 
it  should  be  pointed  out  that  the  need  for  stiffer  control  becomes  less 
urgent  with  increased  range  length  as  more  time  is  available  for 
recovery  from  initial  errors. 

In  conclusion,  it  is  recognized  that  numerous  variations  and 
ramifications  of  the  concept  presented  are  obviously  possible,  most 
of  which  are  left  to  the  reader's  imagination.  Opportunity  is  taken, 
however,  to  emphasize  the  fact  that  a more  efficient  approach  is 
possible  for  situations  where  space  and  launch  accelerations  will 
permit  superconducting  magnets  on  the  projectile.  In  the  meantime, 
subsequent  to  the  experimental  verifications  promised  above,  it 
appears  that  construction  of  a scaled-down,  developmental  version 
of  the  present  concept  would  be  in  order.  This  could  most  easily 
be  accomplished  in  the  presently  inactive  VKF  Tunnel  J (Ref.  15) 
where  proper  power  is  available  for  the  guideway. 
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Optional  Nose  or  Test  Article 


Figure  4.  Projectile  configuration  for  feasibility  study. 
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Figure  6.  Power,  shielding,  and 
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Figure  7.  Transposed  line  current  segments. 
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b.  Image  currents  for  equivalent  fields  exterior  to  cylinder 
Figure  9.  Magnetic  cylinder  in  presence  of  a line  of  current. 
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Figure  10.  Force  exerted  on  a diamagnetic  cylinder  by  a parallel 
line  current. 
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a.  a/b  = 0.1 

Figure  12.  Radial  positional  forces. 
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b,  a/b  = 0.2 
Figure  12.  Continued. 
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d.  a/b  = 0.5 
Figure  12.  Continued. 
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e.  a/b  = 0.8 
Figure  12.  Concluded. 
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s/b 

a.  a/b  = 0.1 

Figure  13.  Angular  positional  forces. 
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b.  a/b  = 0.2 
Figure  13.  Continued. 
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d.  a/b  = 0.5 
Figure  13.  Continued. 
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Figure  15.  Cylindrical  shell  and  line  current. 
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Figure  16.  Nomenclature  for  determination  of  damping  forces. 
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s/b 

a.  a/b  = 0.1 

Figure  17.  Magnetic  damping  coefficients. 
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c.  a/b  = 0.3 
Figure  17.  Continued. 
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e.  a/b  = 0.8 
Figure  17.  Concluded. 
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Figure  18. 


Magnetic  damping  coefficien 
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Dimensions  (in  meters): 
a 1. 27  x 10" 2 

Bore  6.35  x 10"2 

h 1. 016  x 10” 1 

i 9. 525  x 10"3 

1. 905  x 10" 2 
rb  2. 54  x 10'2 

w See  Text 


Figure  20.  Launch  package. 
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Figure  22.  Computed  flight  path  (6,000  m/sec,  'A  atm). 
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Figure  24.  Computed  flight  path  (4,500  m/sec,  'h  atm). 
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Figure  26.  Computed  flight  path  with  increased  lD  (6,000  m/sec,  % atm). 
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APPENDIX  A 
COMPUTER  PROGRAMS 
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«*  PROGRAM  FOR  B-SEGMENTEO/B-INFINITE 

• * 


M SPECIFIES  THE  NUMBER  OF  SEGMENTS  CONSIDERED  (SHOULD  BE  000) 

May 

SEG  IS  THE  LENGTH  OF  ONE  SEGMENT 

SEG=2.0 
DO  1 Lal »S 
AL"L 

0 IS  THE  DISTANCE  OF  ThE  POINT  OF  INTEREST  FROM  THE  Z-AXIS 

0=SEG/AL 

THIS  DO  LOOP  STEPS  THE  POINT  OF  INTEREST  ALONG  THE  Z-AXlS 

00  1 1*1(41 
E*I  -1 
S=(M-l)/2 

Z IS  THE  Z-COMPONENT  OF  THE  LOCATION  OF  THE  POINT  OF  INTEREST 

Z=(S*.OS*E)*SEG 
RATIO»0,0 
DO  2 N*1«M 
A*N 

C*Z-  A*5EG 
D=Z-(A-1.0>*SEG 

Ro,5*(-l.)*»N*(C/SQRT  <C*#2*B**2)-  O/SQRT  <D**2*B**2>) 

2 RATlOsRATIO  * « 

WRITE (6. 10)M»SEG»BtZ»RATIO 

1 CONTINUE 

10  FORMAT C1H0(3HM  -»I3  ,/lH  « 16HSEGHENT  LENGTH  -tElO.M/lH  »3HB  «♦ 
1E10.A/1H  ,3HZ  =.E10.4/1H  ,24H8-SEGMENTED/B-INF INITE  ».E15.6) 
STOP 
END 
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C • 

C * 

C * 

C * 

c * 
c 
c 
c 
c 
c 
c 

DIMENSION  X(-4liY(4)*R(4),)(C<4).YC(41»A(4>  (4) * ADX (4)  .DERVX (4)  * 

lHL)  Y ( 4 ) * AU  Y 1 4)  « OER  V Y (4 ) »FP  I X (4)  .FP I Y (4)  (DANGLE  ( 100)  .8IGROOO)  , 
2BX  (4) « BY (4 ) (SMRSaD(4) 

PIs3. 1415926535 
OAnGLE<1)3-P1/16. 

HIGR( 1) ■-.05*. 0254 
ARA0=0.2*.l>254 
BRAD* 1 • 0* , 0254 
C 

C (XC(N).YC(NI)  INDICATES  POSITION  OF  NTH  LINE  CURRENT 
C ARAD  (BORE  RADIUS)  t ARAIXMODEL  RADIUS) 

c 

PRINT  25.  RRAD  , ARAD 

25  FURMAT  ( ' 0 BORE  RADIUS  *.  E2U.12.  » MODEL  RADIUS  •.£20.12) 

XC(1)>BRAD 
XC (2) bO.O 
XC(3)b-BHAO 
XC (4) *0,0 
YC( 1 )»0,0 
YC(2)=BRA0 
YC(3)«G.O 
YC (4) a-BR AO 
DO  500  L»2(3 
00  500  Ka2.I4 

8 1GR (K ) *6 IGW (K-l ) ♦, OS*. 0254 
OANGLE  <L ) =OANGLE (L-l ) .PI/16, 

C 

C S (BETA)  DEFINES  THE  MOOEL  POSITION 
C (X(N).Y(N))  INDIGATES  POSITION  OF  NTH  IMAGE  CURRENT 

c 

S»BIGR(tO 
BETA»DANGLE (L) 

0»S**2-2,*S*BRAD*COS(8ETA) +BRA0**2 
X < 1 ) »S«COS (BETA) ♦ ( APA0**2« (BRAO-S«COS  (BETA)) 1 /B 
Y(1)«S«SIN(BETA)-(ARAD**2»S*SIN(BETA) >/B 
CHS**2*6RA0**2-2.«S*6RA0*SIN (BETA) 

X (2) »S*COS(BETA) * (ARAQ**2*S*C0S(8ETA) ) /C 
C 

c 

c (PROGRAM  CONTINUED  ON  NEXT  PAGE) 
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Y(2)«S*SlN(bETA> ♦<AR4D»«2*IBRAD-S*S1N<HETA) ) )/C 
D*5**2»HRaO**2+2 • *S*bR A0*CUS  (BETA) 

X ( 3) »S*C0S  (BET A ) - ( ARA0**2<*  (bRAL)*S*COS  (BETA)  > > /D 

Y (3>sS*SlN(8ETA>-(AHAU**2*S*SIN(BETA)  ) /D 

E=5**2*dRA0**2*2.*S*«RA0*SIN<dETAj 

X (4)  *S*COS (BET A > - < ARA0**2«S*COS (BETA) ( /E 

Y (4) »S*SI W (BET A > “ ( ARA0**2* (bRA0*S*5IN(BETA) > > /£ 

M IS  THE  NUMBER  OF  LINE  CURRENTS  PRESENT  (E.O.  A IN  THE  CASE  OF  A 
QUADRUPOLE  FIELD) 

M=>4 

00  10  JsltM 

00  1 J=1*M 

SMRSQD(I>=*(XC(J>-X(I  ) l»»2MYC(J)-Y(I)  ) **2 

dx<i)»(-i.)4«a-n*(YC(j)-uin/SMRSQD(i)  *2.e-t 

1 BY(I)  = <-1.)**(T-1)*(XC<J)-X(I)  l/SHRSODd)  *2.E-T 
1F(m-1  ) 200  *200*300 

BX(N) — X COMPONENT  OF  H PER  I 
8 Y < N) — Y COMPONENT  OF  H PER  I 

00  DO  2 J=2.M 

BX(I)=0X(I)  *BXU-1  ) 

2 BY(I)»BYU>*BYU-1> 

00  CONTINUE 

FPIX(J>=  BY(M)  *<-!.)**( J*l) 

ID  FP I Y ( J ) *8X 1 M > *(-!.)**( J + l) 

IF  t M— 1 ) 51,51*50 

FP1XT  IS  THE  X COMPONENT  OF  FORCE  PER  I SQUARED  FOR  A 1 METER  LONG  MODEL 
FPIYT  IS  THE  Y COMPONENT  OF  FORCE  PER  I SQUARED  FOR  A 1 METER  LONG  MODEL 

51  FPIXT— FPIX(l) 

FPIYT—FPIY  (1 1 
GO  TO  52 

SO  FPIXT— FPIX  (1>-FPIX(2>-FPIX  (3)  -FPIX<4) 

FPIYT—FPIY  (1}-FPIY(2)-FPIY{3)-FPIY<4) 

52  CONTINUE 

FR*FP1 XT*COS (BETA) *FPI YT*SIN (BETA) 

FBETA— FPIXT*SIN(BETA)*FPIYT*COS(bETA) 

PRINT  100,  S»  BETA,  FPIXT,  FPIYT,  FR,  FBETA 
100  FORMAT ( 1HO , 1HR, l X»E 1 Q.4  *2X  »4HBE  TA»1X»E10,4,2X»2HFX*E10»4,2X, 

22HFY ,E10.A,2X,2HFR  * 1X»E1 0 , A, 2X , 5HF8ETA , 1X,£10.A) 

500  CONTINUE 
STOP 
END 
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* 

* 


PROGRAM  FOR  DAMPING  CONSTANTS 


DIMENSION  JJONT  (4)  .OONUT  (4)  » ANGLE  (A)  .SMLR  <4)  t THETA <4)  « BANGLE (4)  • 

2 A (4*365)  .0<<*,365)  i C<4)  ,0(4)  ,E <4) ,EC<4> .F<4> .G (4,365) » 

3H 14,  365 ) »GT (365) .Hi ( J65) ,SQIR<365)  ,dlGR{99) ,P«RH(365) , PWRB(365> * 
4S'11H(36S>  ,DANGLE(40)  .PHI  (365)  ,U0T(4)  ,SSRSQi)<4> 

P I =3* 1415926535 

ARAO  IS  THE  MODEL  RADIUS,  BRAD  IS  THE  GUIDEmAY  BORE  RADIUS 
S(=BIGH)  AND  BETA(=UANGLE)  DEFINES  The  MODEL  POSITION 
M INDICATES  THE  NUMBER  OF  LINE  CURRENTS  PRESENT 

8RA0»3, 5355*. 0254 

ARAU= ( 0 .5- ( . 049/2. ) ) *.0254 

PRINT  35*  ARAD. BRAD 

35  FORMAT  CO  MODEL  RADIUS  •»  £20.12.  ♦ BORE  RADIUS  *.  E20.12) 

BIGW  < 1 ) a- (BRAD/ 1 0 * ) 

DANGLE(1)=-PI/16. 

DO  500  J=2.S 

DO  500  K=2,10 

01 OR (K)mBIGR(K-I) —6 I GR ( 1 ) 

DANGLE <U)=0ANGLE( J-l) *Pl/lb. 

S=8IGR(K) 

BETA=OANGLE( J) 

M»4 

SIMPSON  INTEGSATI 0N--N  BEING  THE  NUMBER  OF  INTERVALS 

00  1 N-I.360 
IF  (N-ll  50.50.51 

50  PHI (1)30. 

GO  TO  52 

51  PHI  <N)oPHHN-l>*2.*PI/360. 

52  CONTINUE 
1 CONTINUE 

00  55  N=»l  .360 
00  2 I-l.M 
00T ( I ) =1-1 
OOnT ( I ) "DOT ( I ) /2 . 


(PROGRAM  CONTINUED  ON  NEXT  PAGE) 
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DONUT (II ®OONT( I )*PL 
ANGLE  ( I I »[10NUT  ( 1 > -BETA 

SMLR  1 1 ) *SQk T (RRaO**2  ♦S*»2-2.*S*BRAD*C0S<  ANGLE < I > > I 
THETA  ( I)  =OONUHI)  ♦A5I''I(S/SMLH<I)*SINUNGLEU)  > > 

GANGLEIIla  (PHI  <N)  “THET  A < i ) ) 

S5R50D ( I ) »APAO**2  + SMLH ( I > *#2  -2.*ARAO*$MLR  < I > »COS (DANGLE ( I > > 
A(IiN)*2.E-T*(ARA0*SMLW(I  >*COS (BANGLE < I)  >-ARA0«*2>/ 

2 ( SMLR ( I ) *SSRSQD { I > ) 

B(XiN)=  -2.E-7*(AHAD*SIN(HANGlE(IM  *SMLR(I>  /SSRSOO<I)> 
C<1>»(S-ORAO*COS( ANGLE (1 ) ) >/$MLH(I> 

o< i>s-brad*sin ( angle < i > > /SMLnm 

EC ( 1 ) ~l , /SORT  < SMLR ( I > **2-S»*2*SlN ( ANGLE (It) **2) 

E(I)=EC(I ) * (S IN ( ANGLt ( I > ) -S*SIN( ANGLE ( I ) ) MS-BRAD* 

2CUS  (ANGLE  ( I ) ) ) /SMLR(I>**2> 

E ( I ) “-EC 1 1 ) • ( COS (ANGLt ( I ) ) -S*HR AO*S I N ( ANGLE < I ) ) **2 

2/SML«( I)**2) 

G( l»N)=A (I »N)*C< I) + B < 1 *N  t *£  < I ) 

2 Htl »N)=A ( I »N)*0(I) ♦ ri ( I *N) *F  < I ) 

C 

C G=CuRRENT  DUE  TO  RADIAL  VELDC I T 7*S l GMA/ INDUCING  CURRENT 
C H*CUR«ENT  OUE  TOANGULAR  YELOC  I T Y#S  IGH  A/ 1 NOUCING  CURRENT 
C 

IEIM-2)  3. <**5 

3 6T(NI*6U«N) 

HT  IN)»H(1»N) 

GO  TO  a 

4 GT  IN)  *G(1*N>-G<2*N) 

HT (N) »H ( l ,N)-H(2.N> 

GO  TO  a 

5 IE (M-A)  6 ♦ 7 ♦ 7 

6 GT(N)a&(I,N)-G(2,N)«G(3»Nt 
HT (N)sH ( 1 . N) -H (2»N) *H (3*N> 

GO  TO  fl 

7 GT (N)aG(l.N) -G (2»N)*G(3»N)“G(A*N) 

HT  (N>=HU.N)-H(2.N)  *H(J*N)  -M(4,N> 

8 CONTINUE 
SUI«(N)=GT(N)**2 

55  SUIB(N)aHT(N)**2 
0"2. ♦PI/360.  »ARAO 
PARR  ( 1 ) *Q/3.*SQIft  ( 1 ) 

PARB ( 1 ) «0/3.*SQI8 ( 1 ) 

PARR ( 360) =Q/3.*SQIH(J60> 

PARS  (360 ) ■Q/3.*S(JIti  ( 360  ► 

00  10  L=2.359 

IF ( (» l ) »*L)  11*11.12 

11  PARH(L)=A.«Q/3,^S0IR(L> 

PARS (L) ■A,«Q/3»*SOI0(L) 

GO  TO  10 

12  PARR(L)a2,»Q/3.*SQIR(L) 

PARB(L)=2.*a/3.*SQia(Lt 

10  CONTINUE 
C 

C KMDR (S (BET A ) ) a PARR 
C 

c 

C (PROGRAM  CONTINUED  ON  NEXT  PAGE) 
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KMDB <S<BETA)  ) *>  PtfHB 

DO  20  La2«  300 
PwRR<L>aPWRR<L-l > +PXRR  <L) 

20  PWRB(L>*PWM«<L-1>+P*RBID 

PRINT  100.  S»  BETA.  P*RR  ( J6U  ) . PWRtJ  { 360 ) 

100  FORMAT  1 1H0.  lHR.lX.ElO.<4*2X»‘*rtHETAt  lX*E10,4.2X.*»rtKMDH«El0.4» 
22X»  fHKMOBETA.ElO.*.) 

500  CONTINUE 
STOP 
END 
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•• 

**  PROGRAM  FOP  GUIDED  FLIGHT 

** 


HEAL  It IO,ttM,KMD,L»LM, LENGTH, M guide 

complex  BC0F,R00T,YT, ALPHA »UCOF , Rl .R2.R3.R4  GUIDE 

DIMENSION  ACOF(a) .flCOFIS) ,CCUF<5> ,C <B ) ,0 ( B) ,ROOT (4) * YCOF (4)  GUIDE 

DIMENSION  DCOF (5) *£COF (b)  GUIOE 

DIMENSION  8N{]1) ,C0F<11)  GUIDE 

DO  200  N = l*4  GUIDE 

ACOF  <N ) -0 . 0 GUIDE 

YCOF(N)=0.0  GUIOE 

8COF(N)»0.0  GUIOE 

HNINIaO.O  GUIOE 

COF<N)aO.O  GUIOE 

200  HOOTCN)»0.0  GUIOE 

DO  201  N=l,5  GUIDE 

201  CCOF<N)=0.0  GUIDE 

YTaO.O  GUIDE 

ALPHAnO.O  GUIDE 

R1»0.0  GUIOE 

R2=0.0  GUIDE 

R3=»0,0  GUIOE 

R4*0,0  GUIDE 

KOUNTaO  GUIDE 

10  KOUNT=KOUNT*1  GUIDE 


VARIABLES  TO  8E  USED 

P (PRESSURE)  tSB  (REFERENCE  AREA)*  L‘(  SEPARATION  BETWEEN  CP  AND  CG)  * 
CMO-CLA-CO0 (AERODYNAMIC  COEFFICIENTS),  RB(BASE  RADIUS),  IO (ROTATIONAL 
INERTIA),  KM-KMO (MAGNETIC  COEFFICIENTS) , LM (SEPARATION  BETWEEN 
CG  AND  CM),  VOUNITIAL  VELOCITY),  YOIINITIAL  ERRROR) , AO ( INITIAL  ANGLE 
OF  ATTACK),  GAMMAOIANGLE  BETWEEN  Z-AXIS  AND  V),  M (MASS  OF  PROJECTILE), 
U (RATE  OF  CHANGE  OF  ANGLE  OF  ATTACK) , INVL (Z-AXIS  POSITION  INTERVAL), 

DT (TIME  INTERVAL  CORRESPONDING  TO  INVL),  LENGTH  (RANGE  LENGTH) 


READ (5, U0)P»SB,L*CLA,CDO»CMQ,RB,IO  GUIDE 
HEAD (5, LlO)KR,KMO,LM* VO, YO » AO .GAMAO » M GUIDE 
READ(S, 120)DT, INVL, LENGTH, U GUIOE 
IF(P.EQ.O.O)  GO  TO  1091  GUIDE 
WR I T£  < 6 « 170 ) GUIDE 


(PROGRAM  CONTINUED  ON  NEXT  PAGE) 
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c 

c 

c 

c 


*R I TE ( 6 1 1 JO ) P . S8  * L » CL A * COO t C«U • RH  * 1 0 
k(HI  IE  <«>*  140  ) KM*KM0«(.M>  VOi  YOi  AOfOAf 

iCNT»l 
IPNT»l 
Y«V0 
1*10 
T“Q.0 
x»o.a 

RHO=Q*0022831»SI5. J7*#P/7hQ.O 
U*RHO*t/##2/2.0 
rOPRsv*GAMAQ 
CU)  ■0«SB*L*  (CLA+COO) 

C (2)  a ( ♦ (CmQ#Q*SR»2.U*h8*HH/V  ) ). 

CIO) 3KM*LM 
C < ^ ) =KMQ*IM 
C <5)»0*S8*CCA 
C<b)aKM 

C ( 7 ) ■ ( • ( Q*S8*CD0/ VI ) 

C(B1«KML) 

D < 1 > = <C (3) *LM*C (4) *VI /I 
0(2)  =V*  <C (2)  + C ( 4)  °LM)  / 1 
Q ( 3 ) =C ( 3) * V/I 
0(4|:V*(CIIKI3)*LH)/I 
0(5) a- MC (6)«LM) / (M«V) ) + ( (C(7)-C(H| ) 

0 (6)  bC  (8) 

0(7)*C(6)/M 
0 ( B) a ( C (5) *C  ( 6 1 *LM ) /H 


GUIDE 

»M*DT tLENOTHfU  GUIDE 

GUIDE 
GUIDE 
GUIDE 
GUIDE 
GUIDE 
GUIOE 
GUIDE 
GUIDE 
GUIDE 
GUIOE 
GUIOE 
GUIOE 
GUIDE 
GUIOE 
GUIDE 
guide 

GUIDE 
GUIDE 
GUIOE 
GUIOE 
GUIOE 

GUIDE 
GUIDE 
GUIDE 


The  ARRAY  CCOF  CONTAINS  The  COEFFICIENTS  of 
equations  for  both  dispersion  and  the  angle 


the  CHARACTERISTIC 
OF  ATTACK* 


c 

c 

c 

c 

c 

c 


CC0F(1)»V*D(6) 

CCOF  <2>  3-D (2) ♦D<5)*V*0<B) 

CC0F(3)30c<n-0(5)*D(2)-D(7)*V-0(b)*Dm 

CCOF  (4)  sj  (5)  »r)  (4)  »D  ( 7 > *0  ( 2)  ~0  ( 6)  *0  (31  ~0  (8)  *U  4 1 ) 

CCOF (51  3-D < 8) *D(3>-D (71*0(41 

DO  20  Ksl,5 

20  DCOF (K) "CMPLX (CCOF (Ki /CCOF  ( 1 1 »l>.  U> 

CJLL  ?fN0C(I.0C0r«,  lUCOF  (31  ,OCOF  (4)  »DC0F(5>  »R1  tR2»R3*R4> 

ROOT (21 3R2 
ROOT (31 »R3 
ROOT (41 »R4 

203  WRlTE(6.2021Rl»R2tR3,H4 

202  FOWMAT('OR1=  » . 2E 1 2. 3 . ' 0R2*  «.2E12.3,*0  R3=  »»2E12.3.»0  R4- 

lc i 31 

204  CONTINUE 

*RITE<6.I501 

YODPRaD (T)*Y0*Q(8)*A0-D(51 *Y0PR*0 (6) *U 


GUIOE 

GUIOE 

GUIDE 

GUIDE 

GUIDE 

GUIDE 

GUIOE 

GUIDE 

GUIDE 

guide 

GUIOE 

GUIOE 

GUIOE 

♦ * guide 

GUIDE 

GUIDE 

GUIDE 

GUIDE 


THE  ARRAY  YCOF  CONTAINS  THE  COEFFICIENTS  OF  THE  NUMERATOR  OF  THE 
LAPLACE  EQUATION  FOR  DISPERSION. 

YCOF (I)3(v«0(6) !*YQ 


(PROGRAM  CONTINUED 


ON  NEXT  PAGE  1 
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C 

c 

c 

c 


YCOF  <2>a<-0<2)*U<8)*V*0{5)  )*YQ*(V+D(6)  ) *YOPR  GUIDE 

YCOF ( 3>«<0(4)-D(2)*D(5>-D(6)*U(l) t *Y 0 ♦IK 8) *V*AU ♦ I «D ( 2>  *0  < 8) ) *YOPR ♦GUIDE 
•0 (6) *YOOPR+D ( 6) * V*U  GUIOE 

YCOF (4)  a ( Q <4) *0  < 5) -0 (8 ) *D ( 1 ) )*Y0  + (-  D ( 8) *D (2) -D (4 ) *D (61 ) *A0*O (4) *YGUIOE 
*0PR*O<8) *YO0PR*D (8) *V*U  GUIOE 

THE  ARRAY  ACOE  CONTAINS  THE  COEFFICIENTS  OF  T HF  NUMERATOR  OF  THE 
LAPLACF  EQUATION  FOR  THE  ANGLE  OF  ATTACK. 

ACOF(1)»IV*0<6))*AO  GUIOE 

ACOF  (2>*(0<5>*V-0(2)  )*A0-O( O*Y0+U(S)*YOPK*YUUPH^V*U  GUIOE 

ACOF  <3)  = l-D<S)*U(2>-UU)*D(6>-U(7)#V)  *AO*U  (3)  *YU+OU  > *YQPR+0  (5)  *YOGUlDE 
*OPR-0 ( 7) »YOP«*D(S)*V*U  GUIOE 

ACOF  (4>=<-<-D(3)*D(5)-D<l)*D<7) )*Y0-(-D(2)*D(  7>*D(3)*D(6)  ) *A0 ♦DO) GUIDE 
#*YOPR-0( O*Y0DPR-0<  7)*V»U)  GUIDE 

oO  99  J*1 i 4 GUIDE 

YCOF <J)*YC0F(J)/CC0F(1)  GUIDE 

99  ACOF(  J)aACOF  (J)/TCOFU)  GUIDE 

100  CALL  VERT ( YCOF  « HOOT » T * YT  > GUIOE 

CALL  VERT (ACOF .ROOT »T .ALPHA)  GUIDE 

IF(ICNT.NF.IPNT)  GO  TU  10S  GUIDE 

WRI TE ( 6 » 160 ) T.X.V.YT. ALPHA. GAMAO  . GUIDE 


104  IPNT=IPNT*INVL 
10b  CONTINUE 

I CNT  «icnt ♦ 1 
TaT.UT 

IF(T,GT.0.10>  GO  TO  106 
X=V*T 

IF(A-LENGTH) 100. 100.107 

106  WRITE  16.190) 

GO  TO  108 

107  WRITE (6.190) 

108  CONTINUE 
GO  TO  10 

1091  CONTINUE 
110  F ORma  T (8E 10  »3> 

120  FORMAT (£10, 3.13. 2F 10.3) 
130  FORMAT  1 1H0.9HP 

* METERS/lH  .9HL 

* 1 / RADIANS/1 1H  . 9HC00 


GUIOE 
GUIDE 
GUIDE 
GUIDE 
GUIOE 
GUIDE 
GUIDE 
GUIOE 
GUIDE 
GUIDE 
GUIOE 
GUIOE 
GUIDE 
GUIDE 
GUIDE 

.E10.4.GH  MM  HG/1H  ,9HS8  * .E10.4.I0H  SQGUIOE 

= .E10.4.7H  METERS/lH  . 9HCLA  ■ .E10.4.10H  GUIDE 
= .E10.4/LH  .9HCMO  = .E10.4/1H  .9HRB  GUIDE 


• * *£10.9, 7H  METERS/lH  . 9HI  = .E10.4.20H  KILOGRAMS-METERS»*2)GUIDE 

140  F ORM AT ( 1 H .9HKM  = .E10.4.14H  NEWTONS/METER/ 1H  . 9HKMD  ■ .ElOGUIDE 

•.4.21H  NEWTON-SECONOS/METEH/1H  ,9HLM  « .E10.4.7H  METERS/lH  i 9HGUIDE 

••VO  = .E10.4.14H  METERS/SECOND/1H  ,9HY0  » .E10.4.7H  METERS/GUIDE 

*1H  .9HALPHA0  o .E10.4.8H  RAUIANS/lM  .9HGAMMA0  a .E 1 0 .4 t ?H«AO I ANS/ IGUI DE 
•H  .9HMASS  * .£10,4. 10H  KlLOGRAMS/lH  ,9H0T  = .E10.4.8H  SECONOGUIDE 

• S/IH  .9HLFNGTH  = .E10.4.7H  METERS/lH  .9HU  =>  .E10.4,l‘5H  RADIANGUIDE 

♦S/SECONO)  GUIDE 

15  0 FORMAT (1H0.5X.4HTIME.9X .8HDIST  ANCE  »7X .6HVELOCI TY.20X » IHY .2SX. 5HALPGUIDE 
*HA . 19X.5HGAMMA/1H  . 4*. 7HSEC0N0S . 8A.6HMETERS * 7x* 1 OHMETERS/SEC. 16X .6GUI0E 
•HMETERS.22X.7HRADIANS. 17X. (HRAOl ANS)  GUIDE 

160  FURMATUH  . (8E15.6)  ) GUIDE 

170  FORMAT  < LH1 . 37* ..8HTUNNEL  G/1H  .29X.24HELECTR0MAGNETIC  GUIOANCE/IHO. GUIOE 


(PROGRAM  CONTINUED  ON  NEXT  PAGE) 
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•24MTHE  INPUT  PARAMETERS  A t^E ) GUIDE 

180  FURMATPO  MAX  TIME  HE ACHEO  *)  GUIDE 

190  FOHMAT ( • 0 MOUEL  HAS  REACHED  TUNNEL  LIMITS  *)  GUIDE 

END  GUIDE 

SUBROUTINE  OANOC {N*BtCtD*E»AltX2*XJ»X4)  QANDC 

IMPLICIT  COMPLEX«8<A-G«u-Z)  QANOC 

COMPLEX *8 1 QANOC 

[= (0 • * 1 • ) QANDC 

IF (N  ,EQ.  0 ) GO  TO  10  QANOC 

A»<  (4,OO*C*£-<0**2*F.)-O**2)/2.O)  QANOC 

Al  = (C*(H*D-4.00«E>/fj»00>  QANOC 

A2=- I <C**3) /27.U0)  QANOC 

A*A*A 1 ♦ A2  QANDC 

8B»  CSQHT  I (A**2>  ♦ < (FiAD-4,0u*E- ( <C**2> /3.00>  ) **J)/2/.00)  QANOC 

A»-A  QANOC 

CALL  CUBRT  ( A «88»R ) QANOC 

PSTARaH*U-4.00«E-C*C/3.UO  QANDC 

RI=-PSTAR/ (3.00*R)  QANOC 

R=CK*fU*IC/J.00>  > QANOC 

P»  CSQRT  < ( 8 4*2/4 , 00 ) -C*R ) QANDC 

PQ=  CSQHT <U»2500*R<**2-L)  QANDC 

A82=.500*8»R-"D  QANDC 

PPQ2*2.00*P*PQ  QANDC 

IF ( CABS  < AB2-PPQ2 ) ,GT.  CABS ( AB2+PPQ2) ) PU»-PQ  QANDC 

PP*<  C AB5 (P ) ► QANDC 

Al=(l. 0*0.0)  QANDC 

B1»(B/2.Q0) *P  QANDC 

Ci*(R/2.Q0)  *-PU  QANDC 

XU<-81*  CSQRT (Bl**2-4.00«Al*Cl> )/(2.00*Al)  QANDC 

X2=(-01-  CSQRTlBl**2-4.00*Al*Cl> )/(2.00*Al>  QANOC 

81« (6/2. 00 1 -P  QANOC 

Cl=  < R/2 . 00  > -PQ  QANOC 

X3»<-81«  CSQRT<01«*2-4.QO«Al*Cl> )/<2.00*Al>  QANOC 

X4=(-B1-  CSQRT<81«*2-4.00*A1*C1) )/<2.00*Al>  QANOC 

RETURN  QANOC 

10  CONTINUE  QANDC 

P=C-(B*«2/3,00>  QANOC 

Q=D-(B*C/3.00)*<  (2.0Q«*B*»3)/27,00>  QANOC 

Zl=-<Q/2.00)* ( CSQRT ( <Q**2/4.00> ♦ IP**3/27,0Q> > > QANOC 

Z2=- ( Q/2.00 ) - ( CSQ«T< (Q**2/4,00>+ IP**3/27,00> > > QANOC 

IF  ( CABS  ( Z 1 1 «GE  • CABS(Z2nZ=Zl  QANOC 

IF  C CABS (22) .GE.  CABS(ZI))Z=Z2  QANOC 

IF ( CABS ( Z ) »EQ.  0.0) Xl*-(B/3.00J  QANOC 

IF(  CABS ( Zl  .EQ.  O.01<2*-<0/3.'OO>  QANDC 

IF  < CABS  < Z I .EQ.  0.0) X3b-(B/3.00>  QANOC 

IF ( CA6SIZI  .EQ.  0.0) RETURN  QANOC 

BBS" (0.00 1 0.00 1 QANOC 

CALL  CUBRT (ZtBBB»Rl ) QANOC 

R*-(P/(3.0Q*R1) ) QANDC 

Bl»-( ,500>* ( (3. 00**. S 1/2.00 )»I  QANOC 

B2«*-< .500) -(  13,00**.S)/2.00)  «I  QANOC 

Xl*-<B/3.00)*Rl+R  QANOC 

X2»-<B/3,00)+B1*R1*W2*R  QANDC 


(PROGRAM  CONTINUED  ON  NEXT  PAGE) 
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XJ=-(B/3.00) +W2«Rl+wl*H  QANDC 

X**=(J  . 0 (J  QANDC 

RETURN  QANDC 

END  QANDC 

SU8R0UTINF  CUHRT < AA »bH * HR ) QANDC 

IMPLICIT  COMPLE  X*8  ( A-G » li-Z ) QANDC 

R€AL*AH»HUtH18,HTH  QANDC 

REAL*4PI»SSSS  QANDC 

COMPLEXES  I QANDC 

990  CONTINUE  QANDC 

QANDC 

11=1  QANDC 

Z1=AA+BB  QANDC 

Z2=AA-S8  QANDC 

IF ( CA8S(72>  .GE.  C AbS  ( Z 1 > > A=Z2  QANDC 

I F ( CABS(Zl)  .GE.  C ArtS ( Z2 ) ) A=Z 1 QANDC 

8=  CONJfiU)  QANDC 

Hi A=  < A*8 ) /2. 00  QANDC 

Hl0=-I*(A-B) /2.00  QANDC 

HTH=  ATAN2(HlB,HlA)  QANDC 

Hs(HU**2  + HlB**2)**.bOO  QANDC 

PI =3. 141592653589  79300  QANDC 

SS5S=3 .00  QANDC 


RR=(H**( 1.00/3,00) ) * ( COS ( (HTH*  (II-l)*2.l)0*PI)/SSSS>  *IM  SIN ( (HTQANQC 
1H* (1 1-1 ) *2.00*PI)/SSSS) ) ) QANDC 

END  QANDC 

SUBROUTINE  VERT(COF,ROOT,TpTF) 

COMPLEX  ROOT  »BCOF  « TFt  BN  GUIOE 

DIMENSION  C0F(4) ,3C0F(4> ,H0UT<4) *8N(4)  GUIDE 

DO  910  1 = 1.4  GUIOE 

910  8N{  I ) = (COE ( 1 ) *ROOT 1 1 ) **3>  + (COE (2) *ROOT 1 1 ) **2) * (COE (3) *ROOT ( I ) ) ♦CGUIOE 


SOE (4) 


GUIOE 


BCOF ( 1 ) *0N ( 1 ) / ( (HOOT ( 1 ) -ROOT ( 2 ) ) * (ROOT ( 1 ) -ROOT (3) ) * (ROOT (1) -ROOT (4GUIDE 


*)  ) ) 


GUIDE 


BCOF  <2>*BN<2) / ( (ROOT (2) -HOOT (1 > ) * (ROOT (2) -ROOT (3) >*<RQ0T (2) -ROOT (4GUIDE 


«)  > ) 


GUIDE 


6C0F (3) =BN (3) / ( ( ROOT ( 3 > -ROOT ( 1 ) ) * (ROOT (3) -ROOT ( 2 ) ) * (ROOT (3) -ROOT (4GUIDE 
•>))  GUIDE 

BCOF (4)=BN(41/(( ROOT (4) -ROOT (1) ) * ( ROOT (4 ) -ROOT (2) ) * (ROOT (4) -ROOT ( 3GUI0E 
•>)>  GUIOE 

TF “6C0F < l ) *CEXP (ROOT ( 1 ) *T ) ♦BCOF (2 ) »CEXP (ROOT (2) *T ) ♦BCOF ( 3) *CEXP  GUIDE 
* (ROOT ( 3 ) *T ) ♦BCOF (4 ) *CEXP ( ROOT (4) *T)  GUIDE 

RETURN  GUIOE 

END  GUIDE 


/* 
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A 

a 


B 

b 

£ 

C 

c 

c 

c 


DO 

La 

ma 


D 

d 


3 

f(s) 

f(t) 

§ 

9 

h 


NOMENCLATURE 

Magnetic  vector  potential,  Wb/m 

Radius  of  conductive  shell  or  projectile  forebody,  m 
Abbreviated  coefficient  (see  Eqs.  (79),  (80),  and  (82)) 
Magnetic  induction,  Wb/m^ 

Distance  from  z axis  to  guideway  line  currents,  m 

Abbreviated  coefficient  (see  Eqs.  (79),  (80),  and  (82)) 

Coefficient  of  circular  harmonic 

Projectile's  aerodynamic  drag  coefficient 

Projectile's  aerodynamic  lift  coefficient 

Projectile's  aerodynamic  damping  coefficient 

Effective  center  of  magnetic  forces 

Projectile's  center  of  mass 

Effective  center  of  aerodynamic  pressure 

Abbreviated  coefficient  (see  Eqs.  (79),  (80),  and  (82» 

Abbreviated  coefficients  (see  Eqs.  (76)  through  (81) 

Abbreviated  coefficient  (see  Eqs.  (79),  (80),  and  (82)) 

Electric  field,  v/m 

Force  per  unit  length,  N/m 

Function  defined  in  Eq.  (66) 

Laplace  transform  at  function  in  terms  of  complex 
quantity 

Function  in  terms  of  time 
Function  defined  in  Eq.  (67) 

Abbreviated  coefficient  (see  Eqs.  (79),  (80),  and  (82)) 
Length  of  conductive  shell,  m 
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I 

A\ 

1 


i 


J 


A 

3 

K 


K 


M 

MD 


K 


Lself 

L . 

mut 

£ 

M 

n 


P 

q 

R 

ji 


A 

a 


sb 


Electric  current  vector,  amp 
Unit  vector  in  x direction 
Area  current  density,  amp/m 
Projectile  moment  of  inertia,  kg-m2 
Unit  vector  in  y direction 
Magnetic  positional  force  constant,  N/m 
Magnetic  damping  force  constant,  N-sec/m 
Function  defined  in  Eqs.  (84)  and  (85) 

Unit  vector  in  z direction 
Self  inductance  of  conductor,  h 
Mutual  inductance  between  conductors,  h 
Moment  arm  for  aerodynamic  forces,  m 
Moment  arm  for  magnetic  forces,  m 
Projectile's  mass,  kg 

Number  corresponding  independently  to  the  nth  circular 
harmonic,  nth  line  current,  n^  current  segment, 
nth  root,  etc. 

Eddy  current  power  per  unit  length  of  shell,  w/  m 
Dynamic  pressure,  N/m2 

Roots  of  characteristic  equation  determined  by 
Eqs.  (79)  or  (80) 

Radial  position  of  projectile's  axis  with  respect  to 
guide  way  axis,  m 

Unit  vector  in  direction  of  $ 

Radial  distance  from  axis  of  projectile  or  other 
reference  to  item  or  point  of  interest,  m 

Radius  of  projectile's  base,  m 

Unit  vector  in  r direction 

Distance  from  line  current  to  item  or  point  of 
interest,  m 

o 

Base  area  of  projectile,  m 
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5 

3 

{L 

t 

V 

w 

ffi 

w 

x 

y 

z 

z' 

or 

0 

(3 

7 

6 

e 

e 

A* 

p 

s 

0> 

0 

u> 


Shielding  ratio  (see  Eq.  (10)) 

Complex  variable 

Guideway  current  segment  length  (see  Fig.  1),  m 
Time,  sec 

Velocity  of  projectile’s  Cg,  m/sec 

Energy  of  eddy  currents  dissipated  per  unit  length 
of  shell,  J/m 

Energy  stored  and/or  dissipated  in  guideway,  J 
Thickness  of  copper  shell  (see  Fig.  20),  m 
Cartesian  coordinate,  m 
Cartesian  coordinate,  m 

Dispersion  of  projectile’s  Cg  from  z axis,  m 

Cartesian  and/or  polar  coordinate  (coincident  with 
the  range  and  guideway  axis),  m 

Offset  z axis  (see  Fig.  7),  m 

Angle  of  attack,  rad 

Angular  position  of  shell  referenced  to  x axis,  rad 

Unit  vector  in  (3  direction 

Angle  between  V and  z axis,  rad 

Phase  angle  of  circular  harmonic,  rad 

Phase  lag  between  inducing  field  and  field  of  eddy 
currents,  rad 

Angular  position  referenced  to  an  x axis,  rad 
Permeability,  h/m 

Volume  resistivity,  ohm-m;  also  density  of  working 
gas,  kg/m3 

Area  resistivity,  ohm 

Potential  of  shell  end  surface,  v 

Angular  position  referenced  to  an  x axis,  rad 

Angular  frequency,  rad/ sec 
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SUBSCRIPTS 

Numericals 

-Independently  identifies  line  currents,  abbreviated 
coefficients  in  Eq.  (81),  etc, 

D 

Damping  field,  currents,  forces,  etc. 

BF 

Basic  guideway  field  ( see  Fig.  1) 

GW 

Guide  way 

i 

Inducing  field  of  guideway  currents  or  other  sources 

M 

Magnetic  (positional) 

MD 

Magnetic  damping 

M S 

Radial  M 

M3 

Angular  M 

MDS 

Radial  MD 

MD3 

Angular  MD 

n 

Independently  refers  to  n*h  circular  harmonic,  n^ 
line  current,  n^1  segment,  n^  root,  etc. 

nf 

Associated  with  image  of  In 

o 

Independently  refers  to  conditions  of  free  space  or 
initial  flight  conditions 

P 

Positional  field,  currents,  forces,  etc. 

9 

Radial 

SM 

Radial  mechanical 

IR  V 

Radial  velocity 

T 

Total 

TF 

Transposed  field  (see  Fig.  3) 

z 

Component  in  z direction 

0 

Angular  component 
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SUPERSCRIPTS 

' Independently  implies  field  of  induced  eddy  currents, 

image  current  for  equivalent  exterior  field,  and  an 
offset  system  of  polar  coordinates,  or  if  centered  over 
variable  implies  first -time  derivative 

" Image  current  for  equivalent  interior  field,  or  if 

centered  over  variable  implies  second-time  derivative 


102 


